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1.0  INTRODUCTION 


Disturbances  in  the  motion  of  reentry  vehicles  during 
boundary  layer  transition  have  been  observed  with  the  assistance 
of  on-board  telemetered  data  for  many  years.  These  disturbances 
are  most  often  observed  as  (1)  a divergence  in  the  angle  of 
attack  history  of  the  vehicle,  (2)  an  increase  in  the  lateral 
rates  and  lateral  accelerations  of  the  vehicle,  and  (3)  a sudden 
decrease,  rather  than  a continued  increase,  in  the  aerodynamic 
frequency  of  the  vehicle.  These  observed  disturbances  have  been 
directly  correlated  with  the  inception  of  boundary  layer 
transition  on  the  heatshield  of  the  vehicle  as  determined  by 
shallow  thermocouple  measurements,  base  pressure  measurements 
and  many  other  transition  monitoring  devices.  This  study  is 
concerned  with: 

1)  Defining  the  source  of  trajectory  disturbances 
which  are  derived  from  boundary  layer  turbulence; 

2)  Assessing  the  ability  of  existing  transducers  to 
measure  these  perturbing  effects; 

3)  Designating  transducer  requirements  if  measure- 
ment deficiencies  exist. 

Analysis  was  applied  to  define  and  rank  the  perturbing  forces 
while  the  transducer  assessment  was  based  upon  a literature 
survey. 
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1.1  Boundary  Layer  Transition  Analysis 


In  an  effort  to  explain  the  source  of  these  reentry  vehicle 
motion  disturbances  resulting  from  boundary  layer  transition, 
calculations  of  the  transitional  hypersonic  boundary  layer  on 
a typical  conical  reentry  vehicle  have  been  performed.  The 
objective  of  these  calculations  was  to  rank  the  relative 
importance  of  various  laminar-transitional-turbulent  boundary 
layer  parameters  by  determining  their  resultant  effect  on  the 
aerodynamic  forces  and  moments  acting  on  the  vehicle.  The 
transitional  boundary  layer  parameters  addressed  in  this  study 
were : 

1)  Boundary  layer  thickening  due  to  transition; 

2)  Increased  surface  blowing  due  to  transition; 

3)  Increased  skin  friction  due  to  transition. 

The  desired  result  of  this  study  was  the  isolation  of  a single 
parameter  associated  with  the  transitional  boundary  layer  on  a 
reentry  vehicle  which  is  the  source  of  the  observed  vehicle 
motion  disturbances.  Another  segment  of  this  study  addressed 
the  question  of  how  this  parameter  can  be  measured  during  the 
reentry  flight  and  the  development  of  instruments  tailored 
toward  this  measurement. 

Attention  was  necessarily  confined  to  a single  reentry 
vehicle  shape  flying  a single  reentry  trajectory.  The  approach 
used  to  arrive  at  a proper  description  of  the  viscous  flow  field 
over  the  vehicle  during  boundary  layer  transition  and  the 
resulting  aerodynamic  forces  and  moments  acting  on  the  vehicle 
was  as  follows: 

1)  Calculate  the  aerodynamic  heating  rate  history  on 
the  vehicle  during  the  reentry. 
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2)  Calculate  the  ablation  response  of  the  vehicle 
heatshield  and  the  surface  temperature,  char  q as 
and  pyrolysis  gas  mass  transfer  rates  into  the 
boundary  layer. 

3)  Calculate  the  viscous  flow  field  over  the  vehicle 
at  an  altitude  representative  of  boundary  layer 
transition  on  the  vehicle.  Consider  realistic 
transitional  front  asymmetries  and  vehicle  angles 
of  attack. 

4)  Determine  the  net  forces  and  moments  acting  on  the 
vehicle  by  integrating  the  pressures  and  shear 
stresses  over  the  vehicle  surface.  Determine  which 
boundary  layer  parameter (s)  are  the  probable  source 
of  the  vehicle  motion  disturbances  occurring  during 
boundary  layer  transition. 

These  studies  are  explained  in  detail  in  Section  2 of  this 
report. 

1.2  Boundary  Layer  Measurements  - A Literature 
Survey 

A review  of  both  open  and  classified  literature  has  been 
undertaken  with  the  purpose  of  accumulating  data  on  techniques 
for  the  measurement  of  boundary  layer  parameters  on  hypersonic 
reentry  vehicles.  The  survey  was  corollary  work  to  the 
studies  outlined  in  Section  1.2  which  examined  the  forces  on 
the  heatshield  during  reentry.  The  results  showed  that  if 
asymmetrical  transition  were  present,  resultant  forces  caused 
trajectory  errors.  Insofar  as  measurements  were  concerned,  two 
parameters  were  especially  noted:  1)  Transition  asymmetry  and 

2)  magnitude  of  gas  pressure  on  the  body  which  directly  caused 
the  body  to  be  displaced  or  to  be  given  an  angular  moment. 
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Having  been  provided  with  these  guidelines,  the  measurement 
literature  survey  was  initiated  by  selection  of  papers  which 
reported  on  experiments,  either  ground  or  flight  tests,  in  which 
boundary  layer  parameters  were  monitored.  It  was  discovered 
that  the  sensors  used  were  as  varied  as  the  boundary  layer 
physical  phenomena  but  that  just  a few  techniques  were  applied 
by  a majority  of  the  investigators.  Table  I is  a listing  of 
boundary  layer  measurands  along  with  the  type  of  sensors  used. 
Pressure,  temperature/heating  rate,  and  hot  wires  were  observed 
to  be  in  prominence  in  the  literature. 

The  severe  environment  of  an  ICBM  reentry  eliminates  the 
use  of  many  of  the  sensors  listed  in  Table  I. 

All  techniques  are  applicable  to  wind  tunnel  tests  as 
noted.  Two  primary  techniques,  hot  wire  anemometers  and  surface 
thin  films  are  not  very  practical  for  flight  rests  due  to  the 
severe  environment.  Thin  films  can  be  used  but  are  limited  in 
altitude  coverage  if  directly  exposed  to  the  flow. 

Attempts  have  been  made  to  use  some  of  the  more  sophisticated 
tehcniques  in  flight  tests;  r.f.  orobes,  interferometers, 
spectrometers,  radiometers,  and  electrostatic  probes  were  flown, 
but  the  equipment  required  tends  to  be  bulky  while  results  were 
not  directly  applicable  to  this  study. 

The  best  and  most  abundant  data  of  transition  phenomena 
were  obtained  with  pressure,  temperature,  and  heating  rate 
transducers;  therefore  this  report  concentrates  principal 
attention  on  these  three  methods. 


TABLE  I BOUNDARY  LAYER  MEASURANDS 


PARAMETER 

MEASUREMENT  TECHNIQUE 

FLIGHT 

TESTS 

Static  Pressure 

Diaphragm  Transducer 

X 

Pressure  Fluctuations 

Diaphragm  Transducer,  Piezo- 
electric Semiconductors 

X 

Surface  Temperature 

* 

Thermocouples,  RTD's 

X 

Heating  Rate  (AT) 

Thermocouples 

X 

Temperature  Profile 

Hot  Wire,  Thin  Films 

Thickness 

Hot  Wire  Anemometer,  Optical 

Mass  Flow  Rate 

Hot  Wire 

Velocity 

Hot  Wire 

Density 

Laser,  Photography 

Electron  Density 

Electrostatic  Probe 

X 

Electrical  Conductivity 

RF  Probes,  Interferometers 

X 

Blowing  Rate 

Recession  Rate  Sensors 

X 

Friction  Coefficients 

Force  Moment  Balances 

Optical  Radiation 

Spectrometer,  Radiometer 

X 

Infrared  Radiation 

Spectrometer,  Radiometer 

X 

* 

Resistive  Temperature  Detectors 


WIND 

TUNNEL 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 
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In  Section  3.0,  detailed  survey  results  are  given  for 
pressure  and  temperature  with  primary  emphasis  on  the  factors 
which  cause  errors.  Section  4.0  contains  conclusions  and 
recommendations  for  future  efforts  in  the  measurement  of 
transition.  Section  5.0  contains  the  references  and  other 
reports  reviewed  which  were  found  to  have  information  pertinent 
to  this  survey. 
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2.0  BOUNDARY  LAYER  ANALYSIS 


2.1  Baseline  Vehicle  Characteristics 

The  vehicle  considered  in  this  study  was  a state-of-the- 
art  type  vehicle  consisting  of  a conical  blunted  cone  with  a 
9-degree  half-cone  angle,  a nose  radius  of  1.65  cm  and  a total 
vehicle  length  of  1.829  meters.  The  total  vehicle  weight  at 
reentry  was  771  kg  and  the  vehicle  possessed  a nominal  static 
margin  of  5%.  The  heatshield  was  carbon  phenolic  on  both  the 
frustum  and  nosetip.  The  baseline  trajectory  consisted  of  a 
10197  kilometer  minimum  energy  trajectory  with  a 91.44  km 
altitude,  reentry  velocity  of  7.163  km/sec  and  a flight  path 
angle  of  22  degrees  below  the  local  horizontal.  Boundary  layer 
transition  was  assumed  to  occur  between  27.43  km  and  21.34  km 
altitudes. 

The  reference  trajectory  used  in  the  present  study  is 

shown  in  Figure  1.  This  trajectory  was  obtained  from  a three- 

degree  of  freedom  point-mass  simulation  of  the  above  vehicle 

reentering  at  the  above  described  initial  conditions.  The 

variation  of  the  aerodynamic  drag  characteristics  of  the  vehicle 

with  free  stream  Mach  number  and  Reynolds  number  was  considered 

for  this  calculation.  The  ballistic  coefficient  at  30.48  km 

2 

altitude  was  98155  N/m  , based  upon  inviscid  drag  only.  High 
altitude  viscous  interaction,  angle  of  attack  effects  and  shape 
change  effects  are  not  included  and  are  believed  to  be  relatively 
insignificant. 

2.2  Aerodynamic  Heating 

The  aerodynamic  heating  on  the  vehicle  was  determined  for 
the  reference  trajectory  of  Figure  1 using  the  ARROW  code.  The 
ARROW  code  was  developed  at  Kaman  and  is  capable  of  calculating 
the  aerodynamic  heating  on  a sphere-cone  at  zero  angle  of  attack. 


Velocity  '-KM/ sec 


FIGURE  1 REFERENCE  TRAJECTORY 


The  code  uses  general  blunt  body  equations  for  both  laminar  and 
turbulent  flow,  accounts  for  nose  bluntness  effects  on  the 
pressure  and  entropy  distribution  along  the  edge  of  the  boundary 
layer,  and  employs  real  gas  thermodynamic  and  transport 
properties  of  air  throughout. 

The  pressure  distribution  and  shock  shape  are  determined  in 
ARROW  through  combined  theoretical  and  empirical  methods.  Over 
the  range  of  conditions  of  interest,  the  shock  wave  shape  can  be 
adequately  approximated  by  the  power  law  curve  of  Gregorek  and 
Korkan  (Reference  1)  in  the  nose  region,  followed  by  a straight 
line  at  the  sharp  cone  shock  angle,  i.e., 

re = rn  a + ^r<2)i'n  f°r  x * xc  (ia» 

where  A = 1.52  + ^5 

n n 

CO 

n = 0.44  + 0.76 


and 

Rg  = Rc  + (x-xc)tan0c  for  x > xQ  (lb) 


Here  0 is  the  sharp  cone  shock  angle  as  computed  from  the 
c 

correlations  of  Simon  and  Walters  (Reference  2)  and  A^  is  the 
bow  shock  standoff  distance,  given  by  Inouye  (Reference  3)  as 


0.78  en 


(2) 
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xr  denotes  the  match  point  between  the  two  fits,  i.e.,  the  point 
where  the  slope  of  the  power  law  fit  is  equal  to  tanGc.  Once 
the  local  shock  wave  angle  is  known,  the  Rankine-Hugoniot 
equations  in  conjunction  with  a normal  shock  table  can  be  used 
to  find  the  entropy  behind  the  shock. 

The  pressure  coefficient  distribution  on  the  nose  (4'  < tj>  ) 
is  obtained  from  a modification  of  Inouye's  work  (op.  cit.) 


d 2 4 

^ = 1 - nsin  <J>  + vsin  <J>  (3) 

pmax 

where  nominally  n = 1.25  and  v = .302.  Back  on  the  conical 
afterbody  the  pressure  coefficient  is  determined  by  a modified 
version  of  the  equation  developed  by  Arne,  et  al  (Reference  4) 


CP  = , .055  _ x— . 1 _ -6. 094 (x-. 48) 2 

Cpc  1+ell (x-1.48)  ’ x2+. 00384 


where 


x . 11 . 


This  equation  reflects  the  under-expansion  on  the  cone  caused 

by  nose  bluntness  effects  followed  by  a recompression  to  the 

sharp  cone  value  C . The  region  between  the  tangent  point  and 

^c 

x = .11  is  then  interpolated  using  a power  law  fit. 

Once  the  pressure  distribution  and  shock  properties  are 
known,  it  remains  to  determine  the  entropy  distribution  at  the 
edge  of  the  boundary  layer.  Because  of  the  curved  bow  shock, 
the  entropy  varies  from  one  streamline  to  another  within  the 
shock  layer,  i.e.,  an  entropy  gradient  exists  due  to  the  presence 
of  vorticity  in  the  flow.  Therefore,  a mass  balance  technique 
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is  used  to  define  the  appropriate  streamline  at  the  edge  of  the 
boundary  layer  for  each  point  on  the  body. 


P 


ooUoo 


2irrp  (5— 6*  ) . 
e e 


An  isentroni.c  extension  is  then  taken  along  this  streamline  to 
the  body  pressure  so  that  all  properties  at  the  edge  of  the 
boundary  layer  can  then  be  computed. 


The  laminar  boundary  layer  equations  are  modified  versions 
of  Lees'  (Reference  5)  local  similarity  analysis 


. 332C  p*{1*  uer 
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Merc  the  * denotes  that  the  properties  are  evaluated  at  Eckert's 
reference  enthaloy  h* 


h*  = . 22h  + . 28h  + .50h 


w 


(6) 


C is  a correction  factor  to  account  for  non-zero  pressure 
gradient  effects  on  the  Dlasius  non-dimensional  enthalpy  gradient 
at  the  wall  (.332  Pr^  ^ ) . Values  of  6 and  6*  as  a function  of  0 
and  other  flow  parameters  are  obtained  from  the  work  of  Walker 
(Reference  6). 

The  turbulent  boundary  layer  equations  are  modified  versions 
of  Walker's  (Reference  7)  particular  solution  of  the  momentum  and 
integral  equations. 


q = 


* *.25  .25 

. 0296C  P 


(Pr*) 


2/3 


(V 


.2 


(h  -h  ) 
r w 


(7a) 


* ..25  2.25+.25H  .25 

p*U*  u„  r 


t = . 0296C 


(SQ) 


.2 


(7b) 


0 = .037 
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p u r 
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(7c) 


where  S 


fS  * *.25  1.25 

p*U*  uor 


ds 


; * *.25  2.25+1.25H  1.25  . 

S =J  p*p*  uQ  r ds 
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h = h + P_ 
r e r 


1/3  Ue 


C = 


.31 <R*) 


.069 


for  R * < 2.52x10' 
e 

for  R^*  > 2.52x10' 


Hero,  the  * denotes  that  the  properties  are  evaluated  at  Eckert's 
reference  enthalpy,  C is  a correction  factor  to  account  for  the 
fact  that  the  Dlasius  flat-plate  skin  friction  values  are  valid 

7 

only  below  a Reynold's  number  of  roughly  10  , and  H is  the 
boundary  layer  form  factor.  These  equations  assume  fully 
developed  turbulent  flow  over  a smooth  surface  with  a virtual 
origin  at  the  stagnation  point.  Values  of  6 and  6*  as  a function 
of  0 and  other  parameters  are  obtained  from  the  work  of  Walker 
and  Schumann  (Reference  8)  . 

The  above  equations  are  solved  in  a step-by-step  fashion, 

starting  at  the  stagnation  point  and  working  back  on  the  body. 

An  iteration  on  the  entropy  is  performed  at  each  body  point 

until  the  mass  balance  is  satisfied.  Thermodynamic  properties 

of  air  are  calculated  using  the  NASA  Ames  real  gas  tables 

(Reference  3).  The  viscosity  of  air  is  based  on  Sutherland's 

equation  for  temperatures  less  than  3333°K,  and  a curve  fit  to 

Hansen's  (Reference  9)  data  at  higher  temperatures.  The  Prandtl 

number  is  taken  as  constant  at  P =0.72. 

r 

Calculations  of  the  aerodynamic  heating  rates  to  the 
reference  vehicle  were  obtained  using  the  ARROW  code  and  the 
reference  trajectory  of  Figure  1.  The  heating  rate  histories 
at  the  stagnation  point,  sphere-cone  tanqency  point  and  rearward- 
most  cone  station  are  shown  in  Figure  2.  These  calculations  are 
based  upon  an  assumed  wall  temperature  of  2778°K.  The  heating 
rates  shown  in  Figure  2 are  laminar  heating  rates  except  for  the 
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rearward-most  cone  station  where  transition  to  turbulent  flow 
at  27.43  km  altitude  is  depicted  with  the  dashed  line.  The 
increased  heating  due  to  turbulent  flow  at  27.43  km  is 
significant,  being  approximately  a sevenfold  increase. 

The  heating  rate  distributions  over  the  body  for  both 
laminar  and  turbulent  flow  at  24.38  km  altitude  are  shown  in 
Figure  3.  Once  again,  these  calculations  are  for  an  assumed 
wall  temperature  of  2778°K  and  the  turbulent  boundary  layer 
has  an  assumed  virtual  origin  at  the  stagnation  point  of  the 
vehicle. 

2.3  Heatshield  Ablation  Response 

As  mentioned  above,  the  ARROW  code  solutions  provided  a 
description  of  the  aerodynamic  heating  environment  of  the 
reference  vehicle  for  the  reference  trajectory  represented  by 
an  entry  velocity  of  7.163  km/sec  and  an  entry  angle  of  -22 
degrees  at  an  altitude  of  91.44  km.  The  ARROW  output  was  also 
utilized  as  input  to  the  REKAP  code  to  determine  the  ablation 
response  of  the  carbon  phenolic  heatshield.  The  REKAP  reaction 
kinetics  ablation  program  is  a one-dimensional  transient 
thermal  conduction  in-depth  pyrolysis  and  surface  ablation 
computer  program.  The  in-depth  pyrolysis  reactions  are  modeled 
by  Arrhenius  rate  expressions  and  the  surface  ablation  can  be 
represented  by  either  an  equilibrium  diffusion  controlled  or 
vaporization- con trolled  reaction. 

When  a typical  inorganic-reinforced  plastic  such  as  carbon 
phenolic  is  subjected  to  atmospheric  reentry  at  hypersonic 
velocity,  a complex  process  known  as  ablation  occurs  which 
sacrifices  the  outer  layer  of  material  while  protecting  the 
bulk  of  the  material  from  thermal  damage.  The  initial  response 
of  the  heatshield  material  to  the  reentry  thermal  pulse  is  the 
decomnosition  of  the  thermosetting  phenolic  resin,  with  the 
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Heating  Rate  —Watts/cm2 


resin  pyrolyzed  to  a mixture  of  light  gases  and  carbonaceous 
materials,  or  char  layer.  As  the  surface  temperature  of  the 
char  layer  increases,  it  begins  to  react  with  oxyqen  in  the 
boundary  layer.  Below  approximately  1000°K,  the  carbon-oxygen 
reaction  is  kinetically  controlled,  being  dependent  solely  upon 
the  forward  reaction  rate  of  carbon  and  oxygen.  Above  approxi- 
mately 1000°K,  the  carbon-oxygen  reaction  is  limited  by  the 
diffusion  of  oxygen  across  the  boundary  layer.  At  very  high 
temperatures,  carbon  undergoes  sublimation/varporization  with 
free  carbon  given  off,  in  addition  to  undergoing  diffusion 
controlled  oxygen  and  cyano  reactions.  The  success  of  the 
ablation  mechanism  is  dependent  upon  forming  high  char  layer 
surface  temperatures  to  reradiate  the  incoming  heat  plus  the 
availability  of  high  heats  of  pyrolysis  in  the  thermosetting 
resin. 

The  REKAP  ablation  program  is  basically  a transient,  one- 
dimensional thermal  conduction  solution  through  a heatshield 
material  containing  a temperature-sensitive  in-depth  reaction. 
The  conservation  of  energy  within  the  heatshield  is  given  by 

PCp  [H  + ''•VT]  = ’•  (KVT)  - Xwjhi  (8) 

where  the  term  w^  represents  the  mass  rate  of  production  of 
species  j by  chemical  reaction  per  unit  volume  per  unit  time,  h_. 
is  the  enthalpy  of  species  j including  chemical,  and  v is  the 
ablation  rate.  This  equation  is  the  normal  diffusion 
equation  with  the  addition  of  the  chemical  source  term  Ewh 
plus  the  term  v*VT  resulting  from  the  use  of  a moving  frame  of 
reference.  The  chemical  source  terms  are  combined  into  a single 
gas  production  term  w^as  and  the  thermal  degradation  process 
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of  the  plastic  can  be  represented  by  a second-order  Arrhenius 
type  rate  controlled  reaction  over  a large  temperature  ranqe 
as  follows 


where,  for  carbon  phenolic 

Pvp  = 1454  kg/m3 
Pc  = 1185  kg/m3 
B^  = 83.4  sec  1 
B2  = l.OllxlO6  sec-1 

= 12058  joules/gm-mole 
E2  - 26796  joules/gm*mole 

An  energy  balance  at  the  surface  of  the  char  layer  provides  the 
heat  flux  boundary  condition  for  the  numerical  solution.  The 
net  conductive  heat  flux  to  the  char  layer  is  given  by 

-X  ft  ' (CH  °eue  • "wn>  <VHw>  • "‘W  ‘ eoTw  (10> 

o 

where  n is  the  transpiration  blockage  coefficient,  assumed  to 
be  0.77  for  laminar  flow  and  0.30  for  turbulent  flow.  The 
surface  emissivity  e for  a carbon  phenolic  char  is  assumed  to 
be  0.85. 

The  mass  transfer  rate  at  the  wall  m , can  be  determined  for 

w 

the  case  of  diffusion  controlled  oxidation  by  employing  the 
relationship  between  heat  and  mass  transfer  in  a laminar 
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boundary  layer  for  unit  Prandtl,  Lewis  and  Schmidt  numbers  as 
proposed  by  Lees  (Reference  10) 


K 


oe 


m 


w 


K peueCH 

ow  o 


+ 1 


(11) 


where  KQe  is  the  mass  fraction  of  oxygen  at  the  edge  of  the 
boundary  layer  (for  air,  K = 0.2314)  and  K is  the  mass 
fraction  of  oxygen  at  the  wall.  For  the  diffusion  limited 
oxidation,  all  of  the  oxygen  at  the  wall  appears  in  the  form  of 
CO  and 


m = 0.1734  p u C„ 
w ^e  e H 


(12) 


As  the  surface  temperature  rises  to  above  approximately 
5000°R,  the  vaporization  rate  of  atomic  and  molecular  carbon 
species,  such  as  C,  C2 , C^,  and  C^,  all  increase  exponentially. 
In  addition,  homogeneous  and  heterogeneous  chemical  reactions 
involving  nitrogen  and  carbon  begin  to  appear  and  form  cyanogen 
(CN) 2 and  the  cyano  radical  CN.  The  REKAP  mass  transfer  rates 
for  a carbon  surface  undergoing  combined  vaporization  and 
diffusion  controlled  oxygen  and  cyano  reactions  are  based  upon 
the  work  of  Scala  and  Gilbert  (Reference  11)  for  a nine 
component  gas  model  using  O,  C>2 , N,  ^ , CO,  CO2 , C,  C^  and  CN, 
with  updating  of  the  transport  and  thermodynamic  properties  by 
Fogaroli  and  Brant  (Reference  12).  Figure  4 shows  the 
variation  of  the  normalized  mass  loss  in  the  diffusion  controlled 
oxidation  and  vaporization  regimes  as  a function  of  surface 
temperature  and  pressure.  The  normalizing  factor  is  the 
diffusion  controlled  oxidation  mass  transfer  rate. 
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FIGURE  4 NORMALIZED  CARBON  MASS  LOSS  IN  AIR 
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Using  the  above-defined  set  of  equations  to  describe  the 
transient  heat  conduction,  in-depth  pyrolysis  reactions  and 
surface  energy  and  mass  balances,  the  REKAP  reaction  kinetics 
ablation  program  solves  the  transient  heat  conduction  equation 
by  stable,  explicit  finite-difference  equations.  The  thermal 
conductivity  and  specific  heat  of  the  carbon  phenolic  ablator 
are  assumed  to  vary  with  temperature  as  shown  in  Figure  5. 

The  ablator  thickness  was  assumed  to  be  1.041-cm  thick  and  was 
integrally  attached  to  a 0.152-cm  thick  aluminum  substructure 
with  a .076-cm  thick  EA  934  adhesive  layer.  The  effects  of 
having  an  underlying  substrate  attached  to  the  carbon  phenolic 
were  virtually  insignificant  since  temperature  rises  in  the 
substrate  at  the  time  of  transition  were  1-2  degrees  K.  As 
mentioned  above,  the  ARROW  heat  transfer  coefficient  and 
recovery  enthalpy  distributions  over  the  reference  vehicle 
during  reentry  were  employed  by  the  REKAP  program  to  determine 
the  heatshield  ablation  response. 

The  results  of  the  REKAP  ablation  analysis  are  shown  in 
Figures  6 and  7 where  the  surface  temperature  and  mass  transfer 
rate  histories  at  three  body  stations  are  shown.  Surface 
temperatures  are  in  the  1100-1400°K  range  during  laminar  flow 
and  jump  to  approximately  3050°K  within  2 seconds  following 
boundary  layer  transition.  Mass  transfer  rates  shown  approxi- 
mately a six-to-eight  fold  increase  due  to  boundary  layer 
transition.  At  the  1.27-meter  axial  station  on  the  body,  the 
mass  transfer  rate  increases  from  .044  to  .29  kg/m  sec  for 
boundary  layer  transition  at  24.38  km.  In  terms  of  the  free 
stream  mass  flux  p^u^,  the  normalized  blowing  rates  at  the 
1.27-meter  station  for  laminar  and  turbulent  flow  would  be 
.00014  and  .00094,  respectively. 
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FIGURE  5 THERMAL  PROPERTIES  OF  CARBON  PHENOLIC  CHAR 


Surface  temperature  '-lO3 


FIGURE  6 REFERENCE  VEHICLE  SURFACE  TEMPERATURE  HISTORIES 
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2.4  Viscous  and  Inviscid  Flow  Fields 


Critical  to  the  work  being  reported  on  here  is  the  postulated 
asymmetric  transition  which  occurs  during  hypersonic  reentry  of 
conical  shaped  vehicles.  The  degree  of  asymmetry  is  not  well 
understood  except  for  some  evidence  derived  from  ground  tests 
which  indicates  a preference  for  transition  to  first  occur  on  the 
leeward  side  of  conical  shapes  at  angle  of  attack.  Additionally, 
asymmetric  surface  roughness  or  manufacturing  irregularities 
may  be  the  source  of  the  postulated  boundary  layer  asymmetry. 

The  objective  of  the  work  contained  in  this  report  is  not  to 
investigate  the  source  or  the  degree  of  asymmetry  which  may  be 
present  during  boundary  layer  transition,  but  instead  to  begin 
with  the  premise  that  an  asymmetry  exists  and  to  attempt  to 
discover  those  parameters  associated  with  transitional  boundary 
layers  which  are  producing  the  excursions  in  the  vehicle  motion 
commonly  observed  during  transition. 

We  v/ill  attempt  to  characterize  the  laminar  and  transitional 
hypersonic  boundary  layer  on  the  reference  vehicle  at  a single 
altitude  of  24.38  km  and  a single  velocity  of  6.870  km/sec.  Both 
non-zero  and  zero  incidence  angles  will  be  considered.  The 
influence  of  wall  mass  transfer  associated  with  the  ablation  of 
the  carbon  phenolic  heatshield  will  be  investigated.  Significant 
limitations  in  the  analytical  techniques  will  be  pointed  out 
where  appropriate  and  recommendations  for  further  study  will  be 
mentioned. 

2.4.1  Zero  Angle  of  Attack  Analysis 

Zero  angle  of  attack  viscous  flow  field  calculations  for  the 
reference  vehicle  at  the  24.38  km  altitude  flight  condition  were 
obtained  using  the  LTBLCEQL  (Laminar  and  Turbulent  Boundary  Layer 
with  Chemical  Equilibrium)  computer  code  (Reference  13) . This 
computer  code  solves  the  compressible  boundary  layer  equations  for  two 
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dimensional  or  axisymmotric  flows  of  either  perfect  gases  or 
reacting  gases  in  chemical  equilibrium.  The  equations  of  motion 
are  cast  in  non-dimensional  form  using  the  Levy-Lees  variables 
and  reduced  to  ordinary  differential  equations.  The  resulting 
parabolic  differential  equations  are  solved  using  finite 
difference  techniques. 

The  governing  equations  for  compressible  boundary  layer 
flows  contain  the  normal  eddy  viscosity  and  eddy  conductivity 
transport  terms  in  both  the  momentum  and  energy  equations , 
respectively.  By  merely  setting  the  eddy  viscosity  and  eddy 
conductivity  terms  equal  to  zero,  nonsimiiar  laminar  boundary 
flows  may  be  computed.  Options  are  available  for  calculating 
either  instantaneous  transition  from  laminar  to  turbulent  flow 
or  transition  over  a prescribed  distance.  The  effects  of  wall 
mass  transfer  are  included  with  injectant  gases  of  either  air, 
helium,  carbon  dioxide  and  argon.  Boundary  conditions  consist 
of  specification  of  the  wall  temperature  and  wall  mass  transfer 
distributions,  the  free  stream  conditions  and  the  boundary  layer 
edge  pressure  and  entropy.  The  actual  boundary  layer  edge 
entropy  is  solved  for  within  the  computer  code  using  an  input 
radial  entropy  distribution  and  streamline  tracking  techniques. 
The  boundary  layer  edge  pressure  is  normally  determired  from  a 
separate  inviscid  flow  field  solution,  as  is  the  radial  entropy 
distribution . 

Defining  the  turbulent  shear  and  the  turbulent  flux  of 
total  enthalpy  in  terms  of  eddy  viscosity  and  eddy  conductivity 
analogies,  the  governing  equations  for  compressible  laminar  or 
turbulent  two-dimensional  or  axisymmetric  boundary  layer  flows 
can  be  written  as 
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where  j = 0 for  two  dimensional  flow 
j = 1 for  axisymmetric  flow 

The  boundary  layer  equations  are  transformed  by  the  Levy- 
Lees  transformation  equations 


dr-  “ We1^  dx 
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With  the  aid  of  the  definitions 
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the  continuity  equation  may  be  written  as: 


where 
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and 


With  the  further  definitions 


F = f' 
q - H/He 

C - f»‘/pe'>e 
e+  = e/v 
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(20) 


the  energy  and  momentum  equations  can  be  written  in  the  following 
standard  form 


cf>"  + A1<J>  * + A2<j>  + A3  + A 4<pf  = 0 


(21) 


where  the  prime  superscripts  denote  differentiation  with  respect 
to  n and  the  £ subscript  denotes  differentiation  with  respect  to 


£. 

For  the  momentum  equation,  <f 
given  by  the  following  equations: 

A = £1  + % _ V_ 

**  i r*  ' A 


= F and  the  coefficients  are 


1 C 


A 


A2  = 6F/AC 


(22) 


A3  = Bpe/PAC 
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A4  = -25F/Ao 


Aq  = C(l+c+) 


A0/c 


For  the 
are  given  by 


energy  equation,  <J)  = g = H/He , 
the  following  equations: 


and  the  coefficients 


A 


1 


V 

A 

o 


(23) 


and  Aq  = Aq/C. 


The  above  equations  are  constrained  by  the  boundary 
conditions  at  the  wall,  n = 0 , and  at  the  outer  edge  of  the 
boundary  layer,  n = ne • The  boundary  conditions  are 


H h 

at  n = 0:  F = 0,  g = ^ = ^;  V = V, 


H 


w 


atn=ne:  F = 1;  g = 1. 
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Two  boundary  layer  transition  models  are  available  in  the 
LTBLCEQL  computer  code.  The  first  transition  model  changes  the 
flow  instantaneously  from  laminar  to  turbulent  flow  at  a 
specified  location  on  the  body.  Heat  transfer  results  obtained 
with  this  model  show  a much  steeper  rise  during  transition  than 
occurs  in  experimental  data.  On  this  basis,  the  instantaneous 
transition  model  is  not  completely  satisfactory. 


The  second  transition  model  provides  for  a regime  of 
transition  between  laminar  and  turbulent  flow.  The  calculation 
is  started  with  the  flow  laminar  and  the  eddy  viscosity  equal 
to  zero.  At  a specified  location  (the  beginning  of  transition) 
the  turbulent  eddy  viscosity  is  gradually  undamped  until  it 
reaches  its  fully  turbulent  value  at  some  prescribed  distance 
downstream  (the  beginning  of  fully  turbulent  flow) . The 
undamping  parameter  used  for  the  eddy  viscosity  is  referred  to 
as  the  transition  intermittency  factor  given  by 


r = 1 - exp (-At2 ) 


where 


A 

K 


0.412 

x-x 


and 


X = 


x,-x 
f o 


(24) 


The  distances  xq  and  x^  refer  to  the  starting  and  ending 
positions  of  boundary  layer  transition. 

For  turbulent  and  transitional  flows,  the  Reynolds  shear 
term  is  given  by  an  eddy  viscosity.  The  eddy  viscosity,  c,  is 
evaluated  using  the  concept  of  a two  layer  eddy  viscosity  model 
consisting  of  an  inner  law,  e^,  valid  near  the  wall  and  an  outer 
law,  e , for  the  remainder  of  the  boundary  layer.  The  inner  eddy 
viscosity  law  is  based  on  Prandtl's  mixing  length  concept  expressed 


as 
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(p'v' ) u'  = pl‘ 


Ia  \2 

I du  i 

\dy]  H i 


where  l is  the  mixing  length.  Prandtl  proposed  that  this  length 
will,  in  its  simplest  form,  be  related  to  the  distance  from  the 
wall,  at  least  in  the  region  of  development  of  turbulence.  His 
proposition  that 

= constant  (26) 

dy 

has  been  tested  under  a variety  of  conditions  and  found  to  be 
quite  adequate  in  the  fully  turbulent  portion  of  the  wall 
region.  Van  Driest  (Reference  14 ) has  provided  a modification 
to  Prandtl 's  mixing  length  hypothesis  to  allow  for  damping  of 
the  Reynolds  shear  stress  near  the  wall.  The  expression  given 
by  Van  Driest  is 


H = Ky  (1  - eY//A) 


A = 26v> 


This  expression  is  modified  further  for  the  case  of  mass  transfer 
at  the  wall,  where  the  damping  term  A becomes 


A = 2 6v 


(exp-v„  / VVH  |/Vxu/P 


The  outer  eddy  viscosity  law  is  based  upon  the  work  of 
Clauser  in  combination  with  the  Klebanoff  intermi ttancy  factor 
and  is  expressed  as 


e = K u 6*  y 
o e l 


I 


where  6*  is  the  incompressible  boundary  layer  displacement 
thickness 


•!■/(' -t) 


dy 


(30) 


and  y is  Klebanoff's  in termi ttancy  factor 


Y 


- [■  + 5'5(s)  ] 


6-.-1 


(31) 


Numerically,  the  inner  eddy  viscosity  increases  with  increasing 
distance  away  from  the  wall  until  it  exceeds  the  outer  eddy 
viscosity.  This  crossover  point  is  monitored  within  the  code 
and  is  used  to  define  the  limits  of  the  inner  and  outer  turbulent 
boundary  layer  regions. 

With  the  boundary  layer  equations  written  in  the  standard 
form  shown  in  Equation  (21) , solutions  are  obtained  using  a 
three-point,  finite  difference  method.  The  solution  procedure 
may  be  fully  implicit  or  fully  explicit.  The  implicit  method  of 
Crank-Nicolson  has  proved  to  be  most  satisfactory.  The  finite 
difference  solution  procedure  has  been  developed  for  a variable 
spacing  of  the  nodal  points  in  the  normal  n coordinate  direction. 
This  permits  a close  spacing  of  points  in  the  region  near  the 
wall  where  the  variation  of  fluid  and  dynamic  properties  is 
greatest . 

The  LTBLCEQL  computer  code  requires  the  specification  of  a 
number  of  boundary  conditions  before  a solution  can  be  obtained. 
As  mentioned  before,  the  REKAP  analysis  supplied  information 
reqarding  the  surface  temperatures  and  pyrolysis  and  char  qas 
mass  transfer  rates  for  laminar  and  turbulent  boundary  layer 
heating.  In  addition  to  this,  boundary  layer  edge  properties 
are  required. 
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The  LTBLCEQL  computer  code  relies  on  the  assumption  that 
the  inviscid  and  viscous  flow  fields  are  uncoupled  and  that 
inviscid  solutions  can  be  used  for  the  specification  of  edge  of 
boundary  layer  pressure  and  entropy  distribution.  This 
assumption  is  valid  as  long  as  the  boundary  layer  displacement 
thickness  is  small  in  comparison  to  the  shock  layer  thickness, 
providing  that  an  entropy  entrainment  type  analysis  is  used. 

We  will  proceed  with  the  approach  that  the  inviscid  and  viscous 
flow  fields  are  uncoupled  and  employ  an  axisymmetric  method  of 
characteristics  solution  to  obtain  the  inviscid  flow  field 
properties.  The  validity  of  the  uncoupled  approach  will  be  re- 
examined upon  completion  of  the  boundary  layer  analysis. 

The  inviscid  pressure  and  entropy  distributions  for  the 
zero  angle  of  attack  viscous  flow  field  cases  were  determined 
using  the  NASA  real  gas  method  of  characteristics  solution 
developed  by  Inouye,  Rakich  and  Lomax  (Reference  15 ) . The  axial 
pressure  distribution  and  radial  entropy  distribution  are  shown 
in  Figures  8 and  9.  Note  the  pressure  over-expansion  and  re- 
compression caused  by  the  nose  bluntness  of  the  body.  The  radial 
entropy  distribution  is  converted  to  a body  surface  distribution 
within  the  LTBLCEQL  code  using  an  iterative  solution  of  the 
entropy  swallowing  by  tracking  streamlines  from  the  shock  cross- 
ing point  to  the  boundary  layer  edge  using  a boundary  layer- 
shock  layer  mass  balance. 

Using  the  boundary  layer  code  conditions  given  by  the  method 
of  characteristics  solution  and  the  wall  temperature  and  wall 
mass  transfer  distributions  given  by  the  REKAP  solutions,  viscous 
flow  fields  were  obtained  with  the  LTBLCEQL  computer  code  for 
four  different  boundary  layer  flows.  They  were: 
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Pressure  --lO5  newtons/m2 


Entropy  S/R 


1)  Fully  laminar  with  wall  mass  transfer 

2)  Fully  laminar  without  wall  mass  transfer 

3)  Transitional  with  wall  mass  transfer 

4)  Transitional  without  wall  mass  transfer. 

For  the  transitional  boundary  layer  flow  calculations,  trans- 
ition was  assumed  to  begin  at  an  axial  length  of  .61  meters  and  end 
at  an  axial  length  of  1.22  meter,  thereby  dividing  the  vehicle 
length  into  a forward  one-third  region  with  laminar  flow,  a 
middle  one-third  region  with  transitional  flow  and  a rearward 
one-third  region  with  fully  developed  turbulent  flow.  The 
selection  of  a transitional  boundary  layer  pattern  described 
above  coupled  with  a fully  laminar  pattern  represented  an 
attempt  to  bound  the  problem  of  asymmetrical  transition  on  the 
vehicle . 

The  wall  mass  transfer  rates  and  wall  temperature  distribu- 
tions were  determined  from  the  REKAP  analysis  described  in 
Section  2.3.  The  choice  of  the  wall  mass  transfer  rates  and  wall 
temperature  distributions  during  transitional  and  turbulent 
boundary  layer  heating  was  somewhat  arbitrary  since  the  actual 
values  are  dependent  on  both  the  boundary  layer  state  and  the 
time  duration  of  the  boundary  layer  state.  The  actual  wall 
temperatures  and  wall  mass  transfer  rates  used  for  both  the 
laminar  and  transitional  LTBLCEQL  solutions  are  shown  in  Fiqures 
10  and  11.  Air  was  chosen  to  simulate  the  injectant  gas.  The 
hiah  surface  temperatures  and  mass  transfer  rates  associated 
with  the  nose  tip  environment  were  not  simulated  in  the 
LTBLCEQL  calculations. 
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FIGURE  10  LTBLCEQL  WALL  TEMPERATURE  DISTRIBUTION 


FIGURE  11  LTBLCEQL  NORMALIZED  WALL  MASS  TRANSFER  DISTRIBUTION 
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Shown  in  Figures  12  through  15  are  the  results  of  the 
LTBLCEQL  fully  laminar  boundary  layer  solutions  with  and  without 
wall  mass  transfer.  Figure  12  depicts  the  boundary  layer  edge 
Mach  number  distribution.  This  figure  is  significant  in  that 
it  shows  the  extent  of  shock  layer  entropy  swallowing  by  the 
boundary  layer.  The  streamline  which  passes  through  a normal 
shock  can  achieve  a maximum  Mach  number  back  on  the  cone  of 
only  3.5,  whereas  a streamline  which  passes  through  a conical 
shock  achieves  a Mach  number  of  12.  These  two  values  represent 
the  extremes  in  the  boundary  layer  edge  Mach  number  on  the  cone. 
The  end  of  body  Mach  number  of  8 in  Figure  12  indicates  that 
nose  tip  effects  are  still  significant  even  through  the  stream- 
wise  distance  is  equivalent  to  113  nose  radii. 

Figure  13  portrays  the  boundary  layer  shear  stress  distribu- 
tion on  the  wall  due  to  skin  friction.  The  effects  of  wall  mass 
transfer  are  much  in  evidence  with  the  transpiration  effect 
reducing  the  skin  friction  by  approximately  20  percent.  Also 
much  in  evidence  is  the  pressure  overexpansion  and  recompression 
aft  of  the  nosetip,  indicating  the  strong  coupling  between  the 
skin  friction  and  the  boundary  layer  edge  pressure. 

The  growth  of  the  boundary  layer  over  the  body  is  shown 
in  Figures  14  and  15  where  both  the  boundary  layer  velocity 
thickness  6 and  the  boundary  layer  displacement  thickness  6*  are 
plotted  as  a function  of  the  streamwise  distance.  The  boundary 
layer  displacement  thickness  6*  is  defined  as 


00 

6*  i (l  - vt)dy 


(32a) 
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Mach  number 


FIGURE  12  BOUNDARY  LAYER  EDGE  MACH  NUMBER, LAMINAR  FLOW 


0.12 


FIGURE  14  BOUNDARY  LAYER  DISPI 
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cm 


or,  in  transformed  variables  as 


and  represents  the  amount  of  displacement  of  the  inviscid  flow 
beyond  the  normal  outline  of  the  actual  body  shape  oroduced  by 
the  presence  of  the  boundary  layer.  The  boundary  layer  velocity 
thickness  5 is  defined  as  the  height  above  the  wall  where  the 
velocity  in  the  boundary  layer  is  0.995  that  of  the  inviscid 
velocity. 

The  inclusion  of  wall  mass  transfer  5b  seen  to  increase 
both  the  displacement  thickness  and  the  velocity  thickness,  with 
the  effects  of  mass  transfer  more  observable  in  the  displacement 
thickness  profile.  The  velocity  thickness  is  rouqhly  a constant 
.15  cm  over  most  of  the  body,  unlike  the  displacement  thick- 
ness which  grows  monotoni cally  with  running  length.  The  relatively 
constant  thickness  boundary  layer  is  the  result  of  the  strong 
entropy  gradient  which  exists  along  the  length  of  the  body.  The 
forward  regions  of  the  body  are  typified  by  hiah  entropy,  low 
density  normal  shock  gas  and  the  entropy  decreases  and  density 
increases  at  the  more  rearward  stations.  This  results  in  an 
increase  in  the  mass  flux  within  the  boundary  layer  with  increas- 
ing running  length  without  a corresponding  increase  in  the 
boundary  layer  thickness. 

The  total  standoff  of  the  shock  wave  from  the  body  at  the 
end  of  body  station  is  4.27  cm  as  determined  from  the  inviscid 
flow  method  of  characteristics  solution.  This  corresponds  to  a 
boundary  layer  displacement  thickness  of  .091  cm  shown  in  Figure 
14.  Thus,  the  boundarv  layer  displacement  thickness  is  seen  to  be 
small  in  comparison  to  the  shock  layer  thickness  for  the  case  of 
fully  laminar  flow,  indicating  that  an  uncoupled  approach  to  the 
viscous  and  inviscid  flow  should  be  satisfactory. 
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Transitional  boundary  layer  calculations  at  zero  anqle  of 
attack  were  accomplished  usinq  the  LTBLCEQL  computer  code 
assuminq  that  boundary  layer  transition  beqan  at  the  .61-meter 
axial  station  and  ended  at  the  1.22-meter  axial  station.  Thus, 
fully  laminar  flow  existed  from  the  nose  rearward  to  the  .61-meter 
axial  station  and  fully  turbulent  flow  existed  from  the  1.22-meter 
to  the  1.83-meter  axial  station  (end  of  body).  Solutions  were 
obtained  for  both  zero  wall  mass  transfer  and  the  turbulent 
wall  mass  transfer  distribution  shown  in  Figure  11.  The 
boundary  layer  transitional  flow  was  modeled  usinq  the  trans- 
ition intermi ttancy  factor  described  above . 

The  results  of  the  LTBLCEQL  transitional  boundary  layer 
calculations  are  plotted  in  Fiqures  16  through  19.  Fiqure  16 
depicts  the  boundary  layer  edge  Mach  number  as  a function  of 
the  streamwise  distance.  This  plot  is  an  effective  measure  of 
the  degree  of  shock  layer  entrouy  swallowing  which  occurs  due  to 
boundary  layer  growth  on  the  body.  As  seen  in  the  fiqure,  the 
boundary  layer  edge  Mach  number  increases  rapidly  during  the 
transitional  flow  region  and  approaches  the  sharp  cone  inviscid 
Mach  number  of  12  at  the  end  of  the  transitional  region.  This 
is  duo  to  the  raoid  increase  in  boundary  layer  thickness  and 
mass  flux  during  the  transition  process.  The  presence  of  wall 
mass  transfer  is  seen  to  have  only  a minimal  effect  on  the  edge 
Mach  number  distribution. 


The  streamwise  distribution  of  the  boundary  layer  wall 
shear  stress  is  pictured  in  Fiqure  17.  The  rapid  increase  in 
wall  shear  stress  is  observed  during  transitional  flow,  with  the 
wall  shear  stress  increasing  by  about  a factor  of  six  during 
transition.  The  double  hump  shape  of  the  shear  stress  distribu- 
tion is  caused  by  the  external  pressure  distribution,  with  the 
first  maxima  a direct  function  of  the  pressure  maxima  at  x = ] 
meter  and  the  second  maxima  corresponding  to  the  onset  of  fully 
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FIGURE  16  BOUNDARY  LAYER  EDGE  MACH  NUMBER 
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Streamwise  distance~meters 


FIGURE  17  WALL  SHEAR  STRESS  DISTRIBUTION, TRANSITIONAL  FLOW 


cm 


FIGURE  18  BOUNDARY  LAYER  DISPLACEMENT  THICKNESS 


developed  turbulent  flow.  The  effect  of  wall  mass  transfer  is 
observed  to  be  almost  of  the  same  importance  as  with  the  ful]y 
laminar  solutions,  with  a 10-20  percent  reduction  in  the  turbulent 
skin  friction  due  to  transpiration  effects. 

The  variation  of  the  boundary  layer  displacement  thickness 
6*  and  the  boundary  layer  velocity  thickness  6 along  the  lenqth 
of  the  body  is  shown  in  Fiqures  18  and  19.  The  displacement 
thickness  begins  to  increase  rapidly  during  the  transitional  flow 
region  and  continues  to  increase  at  about  the  same  rate  throughout 
the  turbulent  boundary  layer  region.  The  turbulent  boundary  layer 
displacement  thickness  is  approximately  an  order  of  magnitude 
greater  than  the  laminar  boundary  layer  displacement  thickness  at 
the  end  of  body  station.  This  can  be  compared  to  an  increase  in 
the  turbulent  boundary  layer  velocity  and  displacement  thickness 
of  about  10  percent  for  the  non-zero  wall  mass  transfer  case 
when  compared  to  the  zero  wall  mass  transfer  case.  Thus  the  most 
important  boundary  layer  thickening  mechanism  appears  to  be  the 
onset  of  turbulence  within  the  boundary  layer  rather  than  the 
onset  of  high  blowing  rates  caused  by  turbulent  heat  transfer. 

The  increased  mass  transfer  rates  associated  with  turbulent 
boundary  layer  ablation  produces  increases  in  the  boundary  layer 
thickness  that  is  an  order  of  magnitude  less  than  that  associated 
with  the  onset  of  purely  unblown  turbulence  in  the  boundary  layer. 

The  turbulent  boundary  layer  displacement  thickness  at  the 
end  of  the  body  of  .91  cm  is  approximately  20  percent  of  the 
end  of  body  shock  layer  thickness  of  4.27  cm.  The  assumption  of 
uncoupled  inviscid  and  viscous  flows  is  being  strained  somewhat 
for  the  situation  of  turbulent  boundary  layer  flow  at  24.38-km 
altitude.  While  the  solutions  presented  here  are  adequate,  the 
24.38-km  altitude  at  transition  most  likely  represents  the  limit- 
ing altitude  at  which  the  uncoupled  approach  can  be  employed. 
Calculations  of  transitional  boundary  layer  flows  at  higher 
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altitudes  should  employ  either  an  iterative  approach  to  the 
inviscid-viscous  flow  field  interface  or  a viscous  shock  layer 
approach . 


2.4.2  Anale  of_  Attack  Analysis 

In  practice,  all  reentry  vehicle  fliqht  tests  possess  non- 
zero anqle  of  attack  values  coincident  with  the  occurrence  of 
boundary  layer  transition  on  the  heatshield.  Coupled  with  this 
is  the  experience  qained  in  ground  tests  which  uniformly  show 
leeward  side  progression  of  transition  on  sharp  and  blunt  cones 
at  anqle  of  attack.  Thus,  anqle  of  attack  effects  may  be 
critical  to  the  understanding  of  the  effects  of  boundary  layer 
transition  on  the  dynamic  behavior  of  reentry  vehicles. 

Included  in  this  effort  was  a study  of  anqle  of  attack  effects 
on  the  reference  vehicle  at  the  reference  altitude  of  24.38  km. 
Angles  of  attack  of  1.5,  3.0  and  4.5  degrees  were  considered. 

The  angle  of  attack  boundary  layer  solutions  were  obtained 
using  the  3-D  boundary  layer  code  (Reference  16) . This  code 
treats  the  full  three-dimensional  steady  boundary  layer  analysis 
of  sharp  and  blunt  cones  at  anqle  of  attack  to  supersonic  and 
hypersonic  flows.  The  program  includes  the  effects  of  surface 
mass  transfer  of  various  foreign  gases  at  the  wall  to  simulate 
ablation  and  also  includes  laminar,  transitional  and  turbulent 
boundary  layer  analysis.  The  turbulent  boundary  layer  has  been 
modeled  using  an  invariant  model  of  three-dimensional  turbulence 
which  employs  the  two-layer  eddy  viscosity  mixing  length  approach. 
An  intermittency  factor  has  been  used  through  the  transition 
regime  to  express  the  probability  of  flow  beinq  turbulent  at  each 
solution  point.  The  governing  boundary  layer  equations  are 
integrated  on  a digital  computer  usinq  a marching  implicit 
finite-difference  scheme. 
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The  throe  dimensional  compressible  boundary  layer  equations 
aro  very  similar  to  those  developed  above  for  axisymmotric  flow 
with  the  added  complexity  of  the  velocity  component  w and  the 
fluctuatinq  velocity  component  w’  acting  in  the  <£  direction. 

The  analysis  employed  in  the  3-D  boundary  layer  code  is  limited 
to  ideal  qas  flows  and  does  not  treat  entropy  gradient  effects. 
The  usual  assumption  of  an  uncoupled  approach  to  the  inviscid 
flow  is  in  effect.  For  completeness,  the  compressible  three- 
dimensional  boundary  layer  equations  are  presented  here  without 
derivation  in  terms  of  mean  physical  variables. 

Continuity  Equation 

f*  (pUr)  + !y  (pVr)  + h (PW)  = 0 (33) 
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Transverse  Momentum  Equation 
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Normal  Momentum  Equation 


^ = 0 
3y 


(35) 


(36) 


Energy  Equation 
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Species  Equation 
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where  V 
is: 


v + p'v'/p.  The  equation  of  state  for  each  species 


(39) 


where  R is  the  universal  gas  constant.  Only  one  species  equation 
is  necessary  since  in  a two-component  mixture  the  mass  fractions 
sum  to  unity: 

Ici  - 1 HO) 

i 

The  viscosity  and  thermal  conductivity  are  related  by  the  Prandtl 
number : 


Pr  = pCp/k 


(41) 


where 


Similarly,  diffusion  and  thermal  conductivity  are  related  by  the 
Lewis  number: 

Le  = p D.f  CD/k  (42) 


The  boundary  conditions  on  the  above  equations  are  as  follows: 
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Momentum^  Equations 
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(43) 


(44) 


(45) 


In  the  derivation  of  the  conservation  equations  the  usual 
assumptions  regarding  the  fluctuating  quantities  have  been 
employed.  These  are: 

1)  The  turbulent  level  is  small  and  therefore  terms 
having  the  mean  square  of  the  velocity  fluctuation 
are  dropped  from  the  equations. 

2)  Molecular  transport  parameters  are  approximated 
by  the  mean  flow  counterparts. 

3)  The  rate  of  change  of  moan  flow  properties  in  the 
normal  direction  is  an  order  of  magnitude  greater 
than  the  rates  of  change  in  the  streamwisc  and 
transverse  directions. 
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The  solution  of  the  governing  equations  requires  the 
expression  of  the  turbulent  shear  terms  and  the  turbulent  flux 
of  total  enthalpy  in  terms  of  the  mean  flow  quantities.  A 
popular  concept  used  to  obtain  these  expressions  is  the  eddy 
viscosity,  eddy-conductivity  analogy  with  the  molecular  viscosity 
and  conductivity  where: 

-pu'v'  = ex  3u/3y 

-pv'w'  = 3w/3y  (46) 

and  -pv'H'  = 3H/3y 

and  where  the  dimensionless  transport  parameters  are: 

Prt  = Cp  c/kt 

Let  " pDtCp/kt 

Windward  Plane  Conservation  Equations 

On  the  windward  plane  of  a cone,  the  transverse  (crossflow) 
velocity,  w,  and  3Pe/9i+>  vanish  due  to  symmetry;  however,  the 
crossflow  velocity  gradient  does  not  vanish  and  still  appears 
in  the  continuity  equation.  Under  these  conditions  the  trans- 
verse momentum  equation  would  vanish  completely  at  the  windward 
plane  where  initial  profiles  are  generated  for  the  remaining 
integration  of  the  governing  equations.  To  avoid  this  problem, 
Moore  (Reference  17)  has  suggested  that  the  transverse  momentum 
equation  first  be  differentiated  with  respect  to  <f>  before 
neglecting  terms  which  vanish  at  the  windward  streamline.  This 
procedure  results  in  the  following  transverse  momentum  equation 
at  the  windward  plane: 
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The  remaining  conservation  equations  are  reduced  to  a quasi-two- 
dimensional  form  at  the  windward  plane  where  w = 0.  The  continuity 
equation  serves  as  the  only  coupling  between  the  transverse 
momentum  equation  and  the  remaining  governing  equations.  For  cones 
at  zero  angle  of  attack  the  transverse  momentum  equation  in  either 
form  vanishes  identically  leaving  a completely  axisymmetric 
problem. 


A more  convenient  form  of  the  governing  equations  for 
numerical  solution  is  obtained  by  introducing  two  stream  functions 
defined  as  follows: 

S'(x,y)  = V57  f (£,n)  (48a) 

<Mx,y)  = V2£/r  g (£,n)  (48b) 

where  £,  n are  the  Levy-Lees  transformed  coordinates  defined  as 
follows : 

x 

£ (x)  =f  pryrurr2dx  (49a) 

o 


n (x, <f>  ,y ) 


(49b) 


This  coordinate  transformation  removes  the  singularity  at  x = 0 
and  stretches  the  normal  coordinate.  Accordingly,  the  transformed 
derivatives  become: 
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3 2 3 . 3n  9 

3x  " Pruryrr  3£  9x  9n 


(50a) 


3_  = 3_  3n  3_ 

3<J>  3<)>  3<J>  3ri 


(50b) 


L = pur/Vn  1- 


(50c) 


Satisfyinq  the  continuity  equation  with  the  above  stream  functions, 
the  following  relations  are  obtained: 


pur  = 


(51a) 


pw  = 


(51b) 


- 3¥  _ 3^ 
3x  3tf' 


(51c) 


Using  Equations  (50a) , (50c)  and  (51c)  results  in  the  following 

expression : 


~2-+rix  6r  f'  + 6q‘  + If  + f + 6 H 

p u \i  r 
r r r 


V + 2£  3f/3£  + f + 5 3q/3<}>  = 0 


V = pvr  V5?/prururr  + nx  6r  f ' + 6g' 


5 = 2C/Prururr  (55) 

Differentiation  of  Equation  (51a)  with  respect  to  y using  Equation 
(50c)  gives  the  expression  for  f': 


f ' 
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Similarly,  differentiation  of  Equation  (51b)  with  respect  to  y 
using  Equation  (50c)  gives  the  exoression  for  g': 


g'  = 


Evaluating  the  momentum  Equations  (34)  and  (35)  at  the  outer  edge 
gives  the  pressure  gradients  as: 


9u„  p w„  9u  p w ~ 

_ e e e e _ e e dr 

peue  9x  r 9$  r 9x 


, 9P  9w  p w 9w  p u w « 

-1  e e ^e  e e *e  e e 9r 

r 9 c})  ~ ^eue  9x  r 9<t>  r 9x 


Using  Equations  (50a)  - (59)  the  governing  conservation  equations 
are  now  expressed  in  terms  of  transformed  variables. 

To  obtain  the  transformed  equations  at  the  windward  stream- 
line, two  new  stream  functions  are  introduced  in  order  to  satisfy 
the  windward  plane  continuity  equation,  as  follows: 

4>  = VJF  f (60 ) 


1|>  = V2Tf  g (61) 

The  transformed  equations  are  then  obtained  in  the  same  manner  as 
for  the  qenoral  three-dimensional  case.  For  a cone  at  zero 
angle  of  attack,  the  system  of  transformed  equations  reduces  to 
a fully  axisymmetric  system  without  a transverse  momentum 
equation . 
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Equations  at  the  Stagnation  Point 


At  the  stagnation  point  of  a blunt  cone,  the  boundary-layer 
equations  have  a removable  singularity.  In  the  limit  as  £ -*■  0 
the  expressions  for  £ and  n are: 

C(x)  = peue  due/dx  x4/4  (62) 

and  y 

D(x,y)  = [2Pe/ue  duQ/dxj  1/2y*  p/pQ  dy  (63) 

o 

Also  at  the  stagnation  point  of  a blunt  cone  the  expression  for 
V in  the  windward  plane  conti^u  ..y  equation  becomes: 
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pv 

du 

eue  dx 
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(64) 


Eddy-Viscosity  Models 

Prandtl's  mixing  length  hypothesis  states  that  the  eddy 
viscosity  is  the  product  of  some  characteristic  length  and  the 
normal  velocity  gradient.  For  two-dimensional  flow  this  concept 
leads  to: 

c = p£*2  | 8ii/3y  | (65) 


Prandtl's  studies  assumed  that  the  eddy  viscosity  should  depend 
only  on  local  eddy  scale  and  on  the  properties  of  turbulence. 

Adams  (Reference  18) extended  this  concept  to  the  three-dimensional 
case  by  assuming  that  the  eddy  viscosity  is  also  independent  of 
coordinate  direction.  In  the  three-dimensional  case  e becomes: 


e 


[(Way)2  + Ow/Dy)2] 


1/2 


(66) 
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which  reduces  to  the  two-dimensional  form  when  w = 0.  This  is 
referred  to  as  the  invariant  turbulence  model . 

The  model  used  in  this  investigation  is  the  common  two- 
layer  inner-outer  model  which  uses  the  Prandtl  mixing  length 
theory  and  the  Van  Driest  or  Reichardt  damping  near  the  wall. 
The  mixing  length  distribution  is  as  follows: 


**  = k*y 


**  = A 


{ 0 < y Xy£/k*} 

< y) 


where 


k*  = 0.435 


A = 0.09 


1/2 

y£  = y when  |^{u2  + w2)/{ue2  + we2  )J  = 0.99 

The  inner  law  is  damped  near  the  wall  so  as  to  yield  the  exact 
laminar  shear  stress  term  at  the  wall.  To  accomplish  this,  two 
different  damping  factors  are  used  in  the  computer  code, 

Van  Driest's  with  local  shear  stress  and  Reichardt' s (Reference 
19)  damping  term. 

Cebeci  (Reference  20)  developed  a mass  transfer  correction 
to  Van  Driest's  inner  eddy-viscosity  law  by  modifying  the 

ic 

damping  constant  A . For  turbulent  flows  with  mass  transfer 
Cebeci  determined  the  damping  constant  to  be 

A = 26  exp  (-5.9  v +)  (68) 


v + = v / (t  /p)1/2 
w w'  w'  1 


I 


Reichardt's  expression  for  the  inner  eddy-viscosity  law 
was  obtained  by  curve  fittinq  experimental  pipe  flow  data. 

The  expression  is: 

ci  = **.  - u-°  tanh  (*!#)]  <6<" 

As  can  be  seen,  this  expression  ddes  not  involve  the  velocity 
gradient  terms.  For  this  reason  it  is  preferred  for  use  in 
numerical  solutions,  since  it  usually  requires  fewer  iterations 
to  converge. 

Following  Equations  (66)  and  (67)  , the  outer  eddy-viscosity 
law  is: 

£o  = *2  y£2  9£/9y 

The  outer  eddy-viscosity  law  is  used  in  conjunction  with  the 
Klebanoff  (Reference  21)  intermittency  factor  assures  a smooth 
approach  of  cq  to  zero  as  y -*•  6.  The  modified  law  is: 

eo  = *2  Yi  Y 9E/9y  (70) 

where  y is  Klebanoff's  intermittency  factor: 

Y = [l  + 5.5  (y/5)6]'1  (71) 


Transition  Models 

Two  models  of  transition  from  laminar  to  turbulent  flow  are 
used  in  the  computer  code.  One  model  is  simply  instantaneous 
transition  to  turbulent  flow,  and  there  really  is  no  transition 
region  or  zone  at  all.  In  the  second  case  a smooth  transition 
to  turbulent  flow  occurs  over  a prescribed  distance.  This 
distance  is  known  as  the  transition  zone  and  is  defined  as  the 
distance  between  the  onset  of  transition  at  x = and  the 
beginning  of  fully  turbulent  flow  at  x = XT  at  some  point  down- 
stream. 
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The  probability  of  turbulent  flow  at  any  point  is  expressed 
by  a model  by  Dahwan  and  Narasimha  (Reference  22)  as: 

I f ( X ) = 1 - exp  [o.412  (2 . 917 ) 2 ( (X-Xt)/(XT-Xt)  ) 2]  (72) 

The  transition  intermittency  factor  is  employed  as  a simple 
multiplier  of  the  eddy  viscosity  in  the  governing  equations  and 
therefore  acts  as  a damping  coefficient  for  the  fully  turbulent 
eddy  viscosity. 

Finite-Difference  Method 


The  finite-difference  method  used  in  our  investigation 
follows  the  method  used  by  McGowan  and  Davis  (Reference  23) 
which  is  similar  to  the  method  developed  by  Dwyer  (Reference  24) 
with  modifications  by  Krause  (Reference  25).  The  method  has 
been  further  modified  to  include  variable  spacing  for  the  normal 
coordinate.  The  accuracy  of  this  method  is  second  order.  The 
method  is  stable  for  negative  transverse  velocities  when  proper 
step  sizes  are  chosen. 

The  momentum,  species  and  energy  equations  are  written  in 
standard  parabolic  form  as: 

2 

A0  « 2 A1  Tn  + A2W  + A3  A4  Ik  A5  94)  “ ° (73) 

9n 

where  w is  the  dependent  variable  in  each  case.  The  coefficients 
Aq  through  A^  are  determined  from  the  transformed  governing 
equations . 

The  derivatives  in  Equation  (73)  are  replaced  with  finite- 
difference  expressions  at  each  of  the  grid  points  in  n , 4'  > £ space. 
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The  finite-difference  form  of  the  Equations  (73)  results 
in  simultaneous  linear  algebraic  equations  of  tridiqonal  form 
which  are  solved  by  a method  developed  by  Richtmyer  (Reference 
26).  The  boundary  conditions  at  both  the  wall  and  the  outer 
edge  must  be  specified  for  this  method. 

Edge  of  the  boundary  layer  conditions  for  the  3-D  boundary 
layer  code  are  established  through  the  usual  software  interface 
with  an  uncoupled  inviscid  flow  solution.  Since  the  3-D  code 
does  not  perform  entropy  gradient  calculations,  only  the 
boundary  layer  edge  pressure  need  by  specified.  The  boundary 
layer  edge  entropy  is  determined  from  a straightforward  normal 
shock-isentropic  expansion  calculation  within  the  code.  The 
boundary  layer  edge  pressure  distribution  may  be  determined  a 
number  of  ways,  from  the  very  simple  modified  Newtonian  or 
tangent-cone  theories  to  the  more  complex  method  of  character- 
istics solutions. 

For  the  cases  considered  here,  real  gas  three  dimensional 
method  of  charac-; eristics  solutions  were  attempted  using  the 
Rakich  (Reference  27)  computer  code.  In  all  cases,  the  method 
of  characteristics  solutions  developed  instabilities  in  the 
pressure  distribution  on  the  leeward  side  of  the  vehicle.  These 
instabilities  were  observed  to  occur  farther  forward  on  the 
vehicle  at  higher  angles  of  attack.  Attempts  at  using  the 
method  of  characteristics  pressure  distributions  with  the 
instabilities  present  as  boundary  layer  edge  conditions  in  the 
3-D  code  resulted  in  instabilities  also  occurring  in  the  boundary 
layer  solutions.  The  boundary  layer  instabilities  were  sufficient 
to  cause  a loss  of  several  meridional  solutions  on  the  leeward 
side  of  the  vehicle,  in  spite  of  a forced  relaxation  of  the 
energy  conservation  criteria.  This  difficulty  led  to  the 
abandonment  of  the  use  of  method  of  characteristics  solutions 
for  specification  of  the  angle  of  attack  pressure  distributions. 


71 


Two  approaches  to  the  specification  of  the  angle  of  attack 
pressure  distributions  remained.  A simple  modifiet  Newtonian 
theory  could  be  used  to  describe  the  surface  pressure  distribu- 
tion on  a blunt  cone  at  angle  of  attack  or  the  analytical 
tables  of  Jones  (Reference  28)  can  be  used  to  describe  the 
surface  pressure  distribution  on  a sharp  cone  at  angle  of  attack. 
With  the  latter  approach  one  would  loose  the  nose  bluntness 
features  insofar  as  they  apply  to  the  pressure  distribution. 

The  most  significant  nose  bluntness  effect,  entropy  layer 
entrainment,  was  lost  from  the  outset  due  to  the  inability  of 
the  three  dimensional  viscous  flow  field  solution  to  properly 
account  for  the  entropy  layer  gradient.  In  actuality,  a better 
approximation  to  the  boundary  layer  edge  entropy  can  be  obtained 
with  an  assumed  sharp  cone  flow  than  with  a blunt  cone  normal 
shock-isentropic  expansion,  since  nearly  all  entropy  swallowing 
was  seen  to  occur  upstream  of  the  transition  reqion  for  the 
zero  angle  of  attack  boundary  layer  solutions.  The  method  of 
Jones  was  therefore  selected  for  the  angle  of  attack  boundary 
layer  edge  pressure  distribution. 

In  the  absence  of  nose  bluntness,  the  pressure  distribution 

is  independent  of  axial  location  on  the  vehicle  and  dependent 

solely  on  the  meridional  location.  The  meridional  distribution 

of  pressure  on  a sharp  9-degree  half  angle  cone  for  angles  of 

attack  of  1.5,  3.0  and  4.5  degrees  are  shown  in  Figure  20. 

Implicit  in  these  results  is  a freestream  Mach  number  of  23.14 

2 

and  a freestream  pressure  of  2801  N/m  . 

Due  to  the  addition  of  a sharp  nose  tip  to  the  basic 
reference  vehicle  for  the  angle  of  attack  viscous  flow  field 
calculations,  alternatives  existed  for  maintaining  either  the 
vehicle  length  or  the  base  diameter  the  same  as  the  blunt  nose 
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FIGURE  20  SHARP  CONE  PRESSURE  DISTRIBUTION  AT  ANGLE  OF  ATTACK 
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zero  anqle  of  attack  calculations.  The  latter  alternative  was 
selected,  thus  increasinq  the  overall  lenqth  of  the  vehicle  to 
1.918  meters  as  opposed  to  the  1.824  meters  blunt  nose  vehicle 
lenqth . 

An  asymmetrical  transition  front  was  employed  for  all  three 
anqle  of  attack  viscous  flow  field  solutions.  The  orientation 
of  the  transition  front  with  respect  to  the  body  was  consistent 
with  the  previous  zero  anqle  of  attack  calculations  and  also 
consistent  with  ground  test  transition  experiments  where  leeward 
side  progression  of  the  transition  front  is  observed.  The 
transition  front  was  a simple  linear  variation  between  complete 
laminar  flow  on  the  windward  side  of  the  vehicle  to  transitional 
flow  on  the  leeward  side  of  the  vehicle.  Leeside  transition 
began  at  an  axial  distance  of  .61  meters  and  ended  at  an  axial 
distance  of  1.22  meters.  Fully  turbulent  flow  existed  aft  of 
the  1.22-meters  axial  station  and  fully  laminar  flow  existed 
forward  of  the  .61-meters  axial  station.  The  transition  front 
was  maintained  constant  at  all  three  angles  of  attack. 

The  wall  temperature  and  mass  transfer  rate  distributions 
developed  previously  in  Section  2 . 3 for  the  zero  angle  of  attack 
viscous  flow  field  solutions  were  applied  directly  to  the  present 
angle  of  attack  calculations.  Specifically,  the  distributions 
for  fully  laminar  flow  in  Section  2.3  were  used  for  the  windward 
side  distributions  and  the  distributions  for  transitional  flow 
in  Section  2.3  were  used  for  the  leeward  side  distributions. 

A linear  variation  with  meridional  angle  was  used  between  these 
two  extremes.  All  angle  of  attack  viscous  flow  field  calcula- 
tions employed  non-zero  mass  transfer  rate  distributions. 

The  results  of  the  anqle  of  attack  three  dimensional 
boundary  layer  calculations  are  shown  in  Figures  21  through  26 
The  first  three  fiqures  show  the  streamwise  variation  in  the 
wall  shear  stress  at  four  different  meridional  angles.  The 
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FIGURE  22  WALL  SHEAR  STRESS  DISTRIBUTION,  o 3.0  DEGREES 
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FIGURE  23  WALL  SHEAR  STRESS  DISTRIBUTION,  « 4.5  DEGREES 
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data  shown  in  these  figures  indicate  that  the  maximum  turbulent 
shear  stress  occurs  along  the  90-degree  meridian,  half-way 
between  the  windward  meridian  at  0 degrees  and  the  leeward 
meridian  at  180  degrees.  Inspection  of  the  raw  data  indicates 
that  the  maximum  is  really  near  ^40-45  degrees  meridionally 
and  is  at  the  very  end  of  the  body.  Inception  of  transition  is 
readily  detected  in  the  shear  stress  distribution,  initially 
occurring  on  the  leeside  and  gradually  moving  around  the  vehicle, 
but  never  reaching,  the  windward  meridian. 

Comparison  of  the  shear  stress  distributions  calculated  for 
zero  angle  of  attack  in  Section  2.4.1  with  those  presented  here 
reveal  two  distinct  differences. 

1.  The  effects  of  the  blunt  nose  overexpansion  and 
recompression  for  the  zero  angle  of  attack 
calculations  are  not  in  evidence  here  due  to  the 
sharp  nose  configuration  employed  for  the  non- 
zero angle  of  attack  calculations. 

2.  Turbulent  boundary  layer  shear  stresses  for  the 
zero  angle  of  attack  calculations  appear  higher 
than  those  calculated  for  the  non-zero  angle  of 
attack  cases.  The  superior  treatment  of  entropy 
gradient  effects  in  the  zero  angle  of  attack 
calculations  is  believed  to  bo  the  cause  of  this 
difference . 

The  streamwise  variation  in  the  boundary  layer  displacement 
thickness  for  the  three  angle  of  attack  calculations  are  shown 
in  Figures  24  through  26.  Once  again,  the  effects  of  the  onset 
of  transition  are  much  in  evidence  with,  in  general,  about  a 
factor  of  five  increase  in  the  rate  of  boundary  layer  growth 
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coincident  with  the  transitional  flow  reqion  along  any  meridian. 

The  boundary  layer  displacement  thickness  does  not  change 
significantly  between  angles  of  attack  of  1.5  and  3.0  degrees  but 
dramatic  increases  are  observed  between  3.0  and  4.5  degrees. 

2.4.3  Sharp  Cone  at  Zero  Angle  of  Attack 

Fundamental  differences  are  seen  in  the  shear  stress  distribu- 
tions obtained  for  the  blunt  cone  zero  angle  of  attack  calcula- 
tions obtained  with  the  LTBLCEQL  computer  code  as  opposed  to  the 
sharp  cone  angle  of  attack  calculations  obtained  with  the  3-D 
boundary  layer  code.  The  blunt  cone  data  do  not  appear  to  match- 
up well  with  the  sharp  cor.e  data  between  0-  and  1.5-degrees 
angle  of  attack.  To  provide  a consistent  set  of  data  for  the 
understanding  of  angle  of  attack  effects,  a sharp  cone  solution 
was  obtained  for  zero  angle  of  attack.  Free  stream  conditions, 
wall  mass  transfer  rates  and  wall  temperatures  were  identical  to 
those  mentioned  above  for  both  the  sharp  cone  at  angle  of  attack 
and  the  blunt  cone  at  zero  angle  of  attack. 

The  results  of  the  sharp  cone  at  zero  angle  of  attack 
calculations  are  shown  in  Figures  27  through  29,  where  distribu- 
tions of  wall  shear  stress,  boundary  layer  displacement  thickness 
and  boundary  layer  velocity  thickness  are  depicted  as  a function 
of  the  streamwise  distance.  Results  for  the  blunt  cone  at  zero 
angle  of  attack  are  shown  on  the  same  figures  as  dashed  lines. 
Differences  between  the  two  solutions  are  very  apparent, 
particularly  with  the  wall  shear  stress  distribution.  No  attempt 
has  been  made  to  try  and  isolate  the  source  of  these  differences. 

It  is  believed,  however,  that  th<.  use  of  the  perfect  gas  equation 
of  state,  constant  inviscid  cone  pressure  and  constant  boundary 
layer  edge  entropy  with  the  sharp  cone  results  versus  the  real 
gas  equilibrium  equation  of  state,  method  of  characteristics 
inviscid  pressure  distribution  and  entropy  swallowing  with  the 
blunt  cone  results  are  the  sources  of  the  differences  in  the  two 
solut ions . 
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FIGURE  27  WALL  SHEAR  STRESS  COMPARISON, 
SHARP  VERSUS  BLUNT  CONE 
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2.5  Aerodynamic  Forces  and  Moments 


A significant  step  in  assessing  the  relative  importance  of 
transitional  boundary  layer  phenomena  on  the  dynamic  behavior 
of  reentry  vehicles  is  to  quantify  the  boundary  layer  parameters 
in  terms  of  actual  forces  and  moments  acting  on  the  vehicle. 

The  boundary  layer  displacement  effect  is  manifested  as  an 
increase  in  the  pressure  on  the  body  which  in  turn  causes  changes 
to  the  axial  and  normal  forces  and  pitchinq  moments  acting  on  the 
body.  The  skin  friction  forces  also  cause  changes  to  the  axial 
forces  and  pitching  moments.  A normal  force  can  be  associated 
with  the  heatshield  ablation  rate  since  the  finite  velocity  of 
the  ablation  products  permeating  the  surface  of  the  body  has  an 
accompanying  thrust.  It  is  the  purpose  of  this  section  to 
quantify  those  effects  related  to  boundary  layer  transition 
which  may  be  a source  of  static  or  dynamic  instabilities  for 
reentry  vehicles. 

The  importance  of  carbon  phenolic  ablation  in  producing 
thicker  boundary  layers  by  itself  has  been  addressed  above. 

The  concern  was  that  the  increased  wall  mass  transfer  rates 
associated  with  turbulent  boundary  layer  heating  may  be  the 
cause  of  significantly  thicker  boundary  layers  and  corresponding 
greater  induced  pressures,  independent  of  the  occurrence  of 
transition  and  the  onset  of  turbulence  within  the  boundary 
layer.  By  zeroing  out  the  mass  transfer  associated  with  heat- 
shield  ablation,  it  was  determined  that  the  onset  of  turbulence 
was  the  principal  contributor  to  increased  boundary  layer 
displacement  thickness  and  increases  in  the  induced  pressure 
following  transition.  The  inclusion  of  wall  mass  transfer  rates 
associated  with  carbon  phenolic  ablation  to  turbulent  heating 
only  added  about  10  percent  to  the  boundary  layer  displacement 
thicknesses . 


' ■* 
* 


1 


i 

' 

'■ 

■ 


I <: 


86 


The  significance  of  ablation  thrust  can  be  addressed  usinq 
a straightforward  approach  to  the  conservation  of  momentum, 
which  is  written  as 


F 


v 


dm 

dt 


where  v is  the  velocity  of  the  ablation  products  at  the  wall  and 

dm/dt  is  the  mass  flow  rate  of  ablation  products  at  the  wall. 

The  mass  flow  rate  for  turbulent  heating  is  discussed  in 

2 

Section  2 . 3 and  reaches  a maximum  of  .176  kg/m  sec  on  the 
rear  of  the  cone.  Assuming  thermodynamic  equilibrium  of  the 
ablation  products  at  the  wall  temperature,  the  velocity  of  the 
ablation  products  is  2.68  m/sec,  giving  a total  ablation  thrust 
of 

F/A  = (.176) (2.68)  = .47  N/m2 

2 

In  comparison  to  the  inviscid  pressure  of  57500  N/m  acting  on 

2 

the  cone,  the  ablation  thrust  of  .47  N/m  is  indeed  small 
and  will  hereafter  be  dismissed. 


The  boundary  layer  displacement  thickness  profiles  will  be 
the  starting  point  for  definition  of  the  boundary  layer-caused 
induced  pressures.  The  induced  pressures  will  be  determined 
using  the  full  tangent-cone  hypothesis,  given  as 


‘ 'P  = 2 y ,y2  _ , . 

P y+1  (Ks  1} 

no  * 


Y(Ks-K)"(  +1) 
Y-l  + 2/K  2 


(74) 


where  K = M^B , B = shock  angle 

K = M (Q  +d6*/dx) , 0 = half-cone  anqle 
s 00  c 

d6*/dx  = boundary  layer  displacement  thickness 
gradient 
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(75) 


I 


Kg  and  K arc  related  by 


's  _ y+1 


(S) 


2 1 


Y+3  k2 


Usinq  the  above  equations,  the  induced  pressures  caused  by 
boundary  layer  thickeninq  can  be  determined.  For  the  blunt  cone 
zero  anqle  of  attack  calculations  for  both  fully  laminar  and 
transitional  boundary  layer  flow,  the  tanqent-cone  method  produces 
the  induced  pressures  shown  in  Fiqure  30.  Keepinq  in  mind  that 


the  rearward  two-thirds  of  the  body  experiences  an  inviscid 

p 

n n r a a * - / 1.1 i ' i a i 


pressure  of  about  57500  N/m  , the  transitional  boundary  layer 
induced  pressure  is  seen  to  be  a siqnificant  fraction  of  this. 


Wo  will  postulate  a hiqhly  asymmetrical  transition  front 
for  the  reference  body  at  zero  anqle  of  attack  in  an  effort  to 
bound  the  problem  to  a certain  extent.  Assume  that  one  meridian 
alonq  the  body  experiences  boundary  layer  transition  beqinninq 
at  the  0. 61-meter  axial  station  and  that  a meridian  180  deqrees 
opposite  to  the  first  has  completely  laminar  flow.  Between 
these  two  extremes  the  boundary  layer  induced  pressure  and  shear 
stress  will  be  defined  accordinq  to  the  followinq  distribution 


AP  = AP  cos 
max 


(Hi 

\ max  / 


At  = T?  + (Tfc  — Tj)  cos 
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FIGURE  30  BOUNDARY  LAYER  DISPLACEMENT  THICKNESS  INDUCED 
PRESSURE  DISTRIBUTIONS 
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x = axial  station  in  meters 

xfc  = transitional  flow  shear  stress  at  x 

x ^ = laminar  flow  shear  stress  at  x 

AP  - difference  between  transitional  and 
max 

laminar  induced  pressures 

The  description  of  <J>  above  produces  the  desired  transition 
front  which  begins  along  a single  meridian  at  x = .61  meters  and 

varies  linearly  with  meridional  angle  until  it  reaches  the  end 
of  the  body  coincidentally  with  the  $ = tt  meridian.  Outside 
of  this  transition  region  laminar  boundary  layer  values 
uniformly  apply. 

Integration  of  the  induced  pressures  and  shear  stresses 
over  the  body  surface  is  performed  to  determine  the  forces  and 
moments  acting  on  the  body.  The  body  center  of  gravity  is 
located  at  a point  61-percent  of  the  length  of  the  vehicle, 
which  is  x =1.10  meters.  Using  the  induced  pressure  and  shear 
stress  distributions  given  in  Equations  76,  the  following 
forces  and  moments  results. 

Normal  force  (6*)  2086  Newtons 

Pitching  moment  (6*)  780  N-m 

Pitching  moment  (i)  212  N-m 

The  pitching  moments  from  both  sources  act  in  concert  and  trim 
the  vehicle  out  at  about  1/2  degree  angle  of  attack.  Furthermore, 
the  transitional  meridian  would  be  pitched  up  to  the  leeward  side, 
which  is  precisely  the  side  favored  for  transition.  Thus  the 
presence  of  any  asymmetry  in  the  transition  front  is  reinforced 
by  the  static  stability  of  the  vehicle.  Additionally,  the 
boundary  layer  displacement  effect  is  predominant  over  the  shear 
stress  effect. 
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For  the  sharp  cone  anqle  of  attack  calculations,  the  integra- 
tion of  the  induced  pressures  and  shear  stresses  over  the  entire 
body  is  a straiqht forward  operation.  The  induced  pressures  are 
once  again  determined  using  the  tanqent-cone  method.  The 
results  of  the  force  and  moment  integrations  are  shown  in  Fiqures 
31,  32,  and  33  as  functions  of  the  angle  of  attack.  Fiqure  31 
portrays  the  normal  force  as  a function  of  anqle  of  attack  which 
results  from  the  boundary  layer  displacement  effect.  This  force 
acts  in  a direction  opposite  to  the  normal  force  associated  with 
the  inviscid  pressure  distribution  on  a cone  at  anqle  of  attack, 
and  is  a consequence  of  the  favored  leeward  side  transition 
that  is  postulated.  This  boundary  layer  induced  normal  force  is 
relatively  constant  until  anqles  of  attack  above  3-degrees  are 
experienced,  where  the  force  begins  to  increase  noticeably.  The 
normal  force  calculated  for  a sharp  cone  at  zero  anqle  of  attack 
is  approximately  one-half  the  normal  force  determined  above  for 
the  zero  angle  of  attack  blunt  cone. 

The  variation  of  pitching  moment  due  to  the  boundary  layer 
displacement  effect  as  a function  of  anqle  of  attack  is  depicted 
in  Figure  32.  Once  again,  this  pitching  moment  is  acting  in 
the  opposite  direction  to  the  inviscid  pitching  moment  due  to 
the  leeside  favored  transition.  The  zero  anqle  of  attack  sharp 
cone  pitchinq  moment  is  about  one-half  that  determined  for  the 
blunt  cone,  and  would  cause  the  vehicle  to  trim  at  a smaller 
anqle  of  attack  than  the  1/2  deqree  mentioned  previously  for 
the  blunt  cone. 

The  shear  stress  pitchinq  moment  is  shown  as  a function  of 
angle  of  attack  in  Figure  33.  At  low  angles  of  attack,  the  shear 
stress  pitching  moment  acts  in  concert  with  the  boundary  layer 
displacement  effect  pitching  moment  and  opposite  to  the  inviscid 
aerodynamic  pitching  moment.  This  direction  sense  is  in  agreement 
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Pitching  moment-^ 


FIGURE  33  WALL  SHEAR  STRESS  PITCHING  MOMENT 
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with  the  blunt  cone  zero  anqle  of  attack  results.  At  about 
1.5-deqrees  anqle  of  attack,  the  shear  stress  moment  is  nulled 
out,  and  above  that  anqle  of  attack  the  moment  takes  the 
opposite  sense  of  direction.  Once  again,  the  displacement 
effect  pitchinq  moment  for  a sharp  cone  at  zero  degrees  anqle 
of  attack  is  approximately  a factor  of  four  greater  than  the 
shear  stress  moment  but  both  are  about  one-half  the  values 
determined  previously  for  the  blunt  cone.  This  latter  differ- 
ence is  believed  to  be  due  to  the  differences  associated  with 
entropy  qradient  effects  between  the  two  types  of  solution. 

The  relative  importance  of  the  boundary  layer  effects  in 
comparison  to  the  inviscid  aerodynamic  forces  and  moments  can 
be  visualized  more  effectively  in  Figures  34  and  35  where  the 
total  forces  and  moments  are  shown  as  a function  of  angle  of 
attack,  with  and  without  the  contributions  due  to  the  sharp 
cone  boundary  layer.  Some  caution  should  be  given  to  these  data 
since  it  has  been  observed  that  the  sharp  cone  boundary  layer 
effects  are  about  one-half  that  seen  for  the  zero  anqle  of 
attack  blunt  cone.  Secondly,  while  the  transition  front 
selected  in  this  analysis  was  chosen  to  represent  an  extreme 
condition  as  far  as  the  geometry  of  the  transition  front  is 
concerned,  it  may  not  represent  an  extreme  as  far  as  the 
important  pitching  moment  is  concerned.  For  instance,  it  is 
apparent  from  inspection  of  the  induced  pressure  distribution 
shown  in  Figure  30  that  a greater  moment  would  result  if  the 
induced  pressure  distribution  for  transitional  flow  were  moved 
rearward  about  .3  meters,  providing  much  qreater  moment  arms 
about  the  x = 1.10-meters  center  of  gravity. 

Inspection  of  Figure  35  indicates  that  asymmetrical 
boundary  layer  transition  effects  on  a sharp  cone  produce  a 
270-400  N*m  decrease  in  the  net  pitchinq  moment  acting  on 
the  vehicle.  The  inviscid  aerodynamic  moment  (which  is 
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FIGURE  34  TOTAL  NORMAL  FORCE  WITH  AND  WITHOUT 
BOUNDARY  LAYER  EFFECTS 
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independent  of  boundary  layer  effects)  is  roughly  1360  N-m  per 
deqree  of  angle  of  attack.  Based  upon  this  comparison  of  the 
pitching  moment  contributions , the  sharp  cone  vehicle  shape 
should  be  stable  above  0.3-degrees  angle  of  attack,  at  which 
point  the  inviscid  aerodynamic  moment  should  exceed  the  destabi- 
lizing boundary  layer  induced  moment.  Once  again,  the  reader  is 
reminded  of  the  conservatism  of  the  sharp  cone  boundary  layer 
induced  moment  of  270-400  N-m.  A more  realistic  maximum 
destabilizing  moment  which  considers  blunt  body  boundary  layer 
flow  and  less  severe  boundary  layer  transition  front  asymmetrical 
geometries  is  probably  in  the  range  of  700-1400  N-m,  which 
raises  the  angle  of  attack  to  achieve  stability  to  0. 5-1.0 
degrees. 

In  conclusion,  it  must  be  emphasized  that  asymmetrical 
boundary  layer  effects  are  important  only  at  small  anqles  of 
attack,  say  less  than  1.0  degrees,  and  cannot  be  significant 
contributors  to  the  in-plane  vehicle  aerodynamic  characteristics 
above  that  point.  This  would  further  suggest  that  if  in-plane 
aerodynamics  are  the  cause  of  the  observed  transient  reentry 
vehicle  excursions  during  boundary  layer  transition,  then  small 
angles  of  attack  are  required  for  the  development  of  these 


excursions. 
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3.0  BOUNDARY  LAYER  MEASUREMENTS  - 


SURVEY  RESULTS 


3 . 1 Pressure  Measurements 

Several  investigators  have  reported  on  experiments  in 
which  measurements  were  made  of  pressure  of  the  flow  field 
boundary  layer  on.  reentry  vehicles  or  wind  tunnel  models.  Some 
the  more  pertinent  data  and  statements  concerning  sources  of 
measurement  inaccuracy  are  summarized  in  this  section. 

The  pressure  exerted  on  a reentry  vehicle  by  the  viscous 
flow  field  exhibits  both  average  (or  static)  and  fluctuating 
(or  acoustic)  properties.  The  measurement  of  these  properties 
are  common  excepting  for  the  frequency  range  over  which  the 
properties  are  measured,  i.e.  the  static  pressure  is  normally 
considered  to  bo  a "steady-state"  or  "average"  value  at  the 
moment  the  measurement  value  is  taken,  whereas  the  fluctuating, 
r.m.s.  (root-mean-square),  or  SPL  (sound  power  level)  is  an 
indication  of  the  rate  and  depth  of  change  of  pressure  in  a 
short  rime  interval  near  the  instant  of  measurement.  Figure  36 

a qualitative  illustration  of  pressure  properties  in  the 
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FIGURE  36  FLOW  CHARACTERISTICS 


Several 

measurements 


factors  contribute  to  the  inaccuracy  of  pressure 

including : 


1)  Transducer  size  and  shape 

2)  Port  aerodynamic  and  erosion  effects 

3)  Temperature  environments 

4)  External  noise  interferences. 


3.1.1  Transducer  Size  and  Shape 

The  results  of  an  investigation  into  the  effects  of  trans- 
ducer size,  shape  and  surface  sensitivity  were  reported  by 
TTnite  (Reference  29)  in  1967  ; he  referenced  earlier  work  by 
'-or cos  after  whom  is  named  the  correction  at  high  frequencies  for 
range  diameter  transducers.  The  Corcos  correction  results  from 
the  random  character  and  short  correlation  lengths  of  pressure 
fluctuations  in  the  boundary  layer.  Figure  37,  extracted  from 
Reference  29,  shows  that  high  frequencies  are  attenuated 
significantly  by  larger  diameter  transducers;  flow  conditions  are 
also  significant  in  this  correction  and  convection  velocity 
should  be  computed  for  each  application.  Shape  and  surface 
sensitivity  are  also  to  be  corrected  for  in  measurements  of  the 
high  frequency  pressure  components.  The  effect  of  non-uniform 
surface  response  (e.g.  an  edge  clamped  diaphragm)  effectively 
reduces  the  transducer  diameter  with  respect  to  the  Corcos 
correction  factor. 

Lewis  and  Dods  (Reference  30)  completed  a-  series  of  wind 
L-unnel  experiments  with  11  flush-mounted  transducers  of  differ- 
ent types  and  manufacturers.  Data  were  obtained  at  three  tunnel 
Mach  numbers  and  plotted  as  shown  in  Figure  38.  Their  comments 
on  the  data  were  that  size  effects  dominated  the  low  frequency 
end  of  the  spectrum,  smaller  diameters  yielding  higher  power 
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FIGURE  38  VARIATIONS  IN' THE  ESTIMATED  POWER-SPECTRAL-DENSITY  MEASUREMENTS 
OF  THE  SURFACE-PRESSURE  FLUCTUATIONS  AS  DETERMINED  FROM  11 
DIFFERENT  TRANSDUCERS  ( REFERENCE  38  ) 
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spocti.nl  density  .levels  at  all  frequencies  than  larger  diameters, 
Transducer  mounting  and  system  noise  levels  were  ruled  out  and 
uncertainties  were  attributed  to  some  undefined  surface  pressure 
Uuctuaf ion-transducer  interaction.  This  conclusion  implies 
potential  errors  may  be  encountered  with  a static  pressure 
measurement,  even  with  restricted  bandwidth. 


1.1.2 


Port:  Ae rodynamics 


A small  diameter  tube  or  port  connecting  the  flow  field 
to  t no  sensor  has  been  applied  by  many  researchers  to  separate 
the  sensor  from  the  boundary  layer  environment.  The  port, 
however,  introduces  some  additional  effects  which  cause  the 
measured  pressure  to  be  in  error.  Heller  and  Widnall  (Reference 
•’  1 ) found  that  the  high  frequency  response  with  flow  is 
increased  over  that  for  no  external  flow,  the  change  depending 
on  flow  field  parameters;  the  data  presented,  however,  was  not 
conclusive  and  was  limited  to  7 kHz  maximum  frequency. 

Johnson  and  Macourek  (Reference  32)  made  measurements  with 
ported,  flush,  and  back-surface  ..transducers  in  a wincLtunnel  at 
M-v>  = 4-  ThGir  data  showed  that  noise  in  the  boundary  layer 
increased  during  transition  and  diminished  immediately  after 
but  continued  to  increase  as  tunnel  pressure  was  increased; 
data  were  taken  up  to  20  kHz.  They  noted  that  acoustic  detection 
of.  transition  is  consistent  with  thermal  methods.  Differences 
between  ported  and  f .lush -mounted  transducer  data  at  low 
frequencies  were  attributed  to  tunnel  noise. 

Hy  using  flush  transducers  for  reference,  Franklin  and 
Wallace  (Reference  33)  measured  the  correction  required  in 
ported  static  measurements  in  wind  tunnels,  using  extreme  care 
to  machine  the  edges  of  the  holes.  They  state  that  even  if  a 
hole  is  near-perfect,  a residual  error  remains  due  to  local 
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changes  in  boundary  conditions  and  consequent  disturbance  of  the 

★ 

boundary  layer.  By  referencing  their  data  to  the  wall  shear 
stress,  they  obtained  a consistent  error  curve  which  also 
included  the  tunnel  flow  conditions;  their  data  are  shown  in 
Figure  40.  It  should  be  noted  that  the  error,  e,  relative  to 
the  shear  stress  is  very  large  but  smaller  in  relationship  to 
the  body  pressure  which  may  be  an  order-of-magn itude  or  more 
greater  than  the  shear  stress. 

Cassanto  has  reported  on  wind  tunnel  and  flight  test 
experiments  using  ported  gages  (References  35,  36,  37,  and  38). 
His  data  showed  that  transition  is  detectable  but  depends  to  a 
certain  extent  on  the  geometry  of  the  transducer  mount.  Figure 
41  illustrates  different  pressure  signatures  obtained  by  differ- 
ent port  geometries  in  a wind  tunnel  environment.  Cassanto 
postulates  that  this  discrepancy  could  possibly  be  attributed  to 
model  misalignment.  Earlier  investigations  (Reference  32)  have 
revealed  that  porting  and  porting  geometry  effects  both  the 
spectral  density  and  the  total  pressure  measured  in  a wind 
tunnel  environment  (Figures  42  and  43).  Cassanto  also 
investigated  port  erosion  effects  in  a rocket  exhaust  facility 
and  showed  that  eroded  ports  record  15  to  30  percent  lower 
pressures  than  non-eroded  ports.  He  also  suggested  a layered 
disk  arrangement  to  reduce  port  erosion  effects  in  an  ablation 
environment . 

Pressure  measurements  on  flight  test  reentry  vehicles  have 
been  partially  successful  usir g ported  gages.  Data  shown  in 
Figure  44  show  ported  sensors  responsed  to  transition  but 

it  ' 

Coe,  in  Reference  34,  reinforced  this  observation  in  a discus- 
sion on  transducer  "flushness"  by  stating  that  very  small 
discontinuities  on  the  order  of  y/5  <.  i.001  significantly 
affected  rms  intensities  as  well  as  coherence  and  convection- 
velocity  measurements  (see  Figure  39)  . 
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( REFERENCE  34) 
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• D IA . = 1 . 83  cm 
A DIA.  = .914  cm 


tf-Flush  tronsducer  data 


Dotted  curves  obtained  with 
transducer  behind  test  static 
holes. 
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a)  PRESSURE  DROP  DATA  OBTAINED  WITH  FLUSH  TRANSDUCERS 


AP*  =Normalized  pressure  difference 
^^Normalized  shear  stress 
e =Pressure  error 
r =Shear  stress 
-Friction  velocity 
d=Hole  diameter 
^“Kinematic  viscosity 
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FIGURE  41  COMPARISON  OF  FLUCTUATING  PRESSURE  DATA  FROM  THE  SHORT 
AND  LONG  PORT  AT  3.3  NOSE  RADII  ( REFERENCE  35  ) 


OVERALL  SOUND  PRESSURE  LEVELS 


GAGE  MEASUREMENT,  FLIGHTIEST  ( REFERENCE  39 ) 


produced  anomalous  results  thereafter  (Reference  39) . In 
another  fliqht  test  (with  a different  reentry  vehicle)  results 
down  to  transition  followed  predictions,  but  were  also 
anomalous  thereafter;  Figure  45  shows  these  data  (Reference  39). 
These  data  also  indicate  discrepancies  which  can  arise  between 
identical  ports  on  the  same  reentry  vehicle. 

Flight  data  obtained  with  different  port  geometries  on 
the  same  flight  show  a marked  difference  in  the  pressure 
signatures  at  onset  of  transition  (Figure  46)  . 

3.1.3  Temperature  Environments 

In  wind  tunnel  tests,  boundary  layer  gas  temperatures 
approach  the  plenum  temperatures;  the  model  will  approach  the 
same  temperature  in  a long  test  period.  Typical  temperatures 
encountered  range  from  480°K  to  1100°K;  thus  sensors  mounted 
flush  with  the  surface  will  encounter  transient  temperature 
effects  which  cannot  be  calibrated  in  the  transducers  presently 
available  from  vendors.  This  source  of  error  was  not  discussed 
in  the  reports  covered  in  the  survey. 

Errors  of  this  type  will  potentially  be  even  larger  during 
flight  tests  where  much  larger  temperature  gradients  are 
encountered. 


3.1.4  External  Noise  Interference 


All  pressure  measurements  in  wind  tunnels  are  made  in  high 
level  noise  environments  due  to  the  noise  generated  by  turbulence 
on  tunnel  walls  or  other  disturbances  upstream  from  the  test 
model.  Reshotko  (Reference  40) , Kendall  (Reference  41)  and 
Beckwith  (Reference  42),  Owen  (Reference  43)  and  Harvey 
(Reference  44)  separately  addressed  the  influence  of  tunnel 
noise  on  transition  in  the  boundary  layers  of  test  models. 
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FIGURE  45  PORTED  GAGE  MEASUREMENTS  (REFERENCE 


There  was  aqreement  that  wind  tunnel  noise  is  of  such  a magnitude 
as  to  be  the  predominant  factor  in  transition  on  test  models. 
Pitot  measurements  by  Harvey  (Reference  44),  et  al  at  Ma=  5 
showed  that  background  fluctuation  levels  reached  as  high  as 
3%  of  p.  with  frequency  components  out  to  250  kHz;  Figure  47  is 
extracted  from  his  report.  Typical  boundary  layer  noise  levels 
ranqe  from  0.2  to  1%  of  p^. 

The  principal  reason  for  citinq  the  above  is  to  illustrate 
that  wind  tunnel  test  data  of  boundary  layer  noise  must  be 
corrected  for  the  high  background  level  or  erroneous  inter- 
pretation may  result.  Also,  in  very  noisy  tunnels  wherein 
pressure  sweeps  are  used  to  facilitate  transition  measurements, 
turbulence  may  be  initiated  erratically  and  provide  apparently 
anomalous  results.  Figure  48  (Reference  42)  illustrates  sources 
of  noise  in  wind  tunnels  and  Figure  49  (Reference  42)  shows  that 
boundary  layer  noise  was  significantly  larger  due  to  turbulence 
only  at  a peak  near  60  kHz  while  tunnel  noise  dominated  below 
30  kHz. 


3.2  Temperature  Measurements 

Classical  definitions  of  transition  evolved  from  tempera- 
ture measurements  of  wall  tempeature  and  surface  heating  rates. 
Using  thermocouples  either  singly  to  obtain  wall  temperature, 
or  differentially  in  pairs  to  obtain  heating  rate,  data  are 
plotted  versus  distance  or  versus  Re^  on  conical  or  inclined 
flat  plates  and  certain  features  are  then  identified  which  are 
related  to  transition.  Recent  plots  by  three  investigators  are 
shown  in  Figures  50,  51,  and  52;  these  data  were  presented  by 
Owen  (Reference  43) , Demetriades  (Reference  45)  and  Martellucci 
(Reference  46)  respectively. 
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FIGURE  47  EFFECT  OF  NOZZLE  WALL  TEMPERATURE  ON  DISTURBANCE 
LEVELS  IN  FREESTREAM  ( REFERENCE  44) 
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a)  VARIATION  OF  OVERALL  RMS  PRESSURE  WITH  UNIT  REYNOLDS  NUMBER 


b)  SPECTRA  OF  PRESSURE  FLUCTUATIONS  WITH  Af  200HZ 


FIGURE  49  MEASUREMENTS  OF  FLUCTUATION  PRESSURES  UNDER  LAMINAR  AND  BOUNDARY 
LAYERS  ON  A SHARP  CONE  IN  MACH  6 HIGH  REYNOLDS  NUMBER  TUNNEL  AT 
NASA  LANGLEY-  ( REFERENCE  42  ) 
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FIGURE  52  HEAT-TRANSFER  DISTRIBUTION  ON  A CONE 
(REFERENCE  46  ) 
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Heatinq  rate  (heat  flux  or  heat  transfer)  data  have  been 
plotted  with  Stanton  number  as  the  dependent  variable  by  all 
three  authors.  The  Stanton  number  is  defined  by  the  following 
relationship: 


noo 


C p V 

P OO  00 


where 


heat  flow  film  coefficient  at  the  wall 
specific  heat 

gas  density 

gas  velocity. 


Transitional  features  are  noted  by  each  author  - onset,  beginning, 
end,  and  peak.  Onset  is  characterized  by  an  upturn  in  the  data 
away  from  the  straight  line  given  by  the  above  equation  where 
pro  and  Vot  .are  increasing.  Beginning  has  been  defined  by  Owen 
to  be  the  point  where  the  laminar  slope  and  transition  slope 
lines  intersect.  The  peak  is  the  maximum  value  of  the  data  and 
the  end  is  the  point  at  which  full  turbulent  slope  intercepts 
the  transition  slope.  It  should  be  noted  that  Martellucci  shows 
"end"  at  the  peak  rather  than  at  the  slope  intercept  (Fiqure  52), 
whereas  Owen,  et  al  show  all  four  points  - a beginning  and  an 
end  which  do  not  correspond  with  either  of  the  other  authors' 
definitions.  (These  comments  are  not  intended  to  be  critical 
of  the  individual  workers,  but  rather  to  illustrate  the  lack  of 
a common  definition  of  transitional  features.) 

Owen  also  shows  (Figure  50)  definitions  of  transitional 
features  based  on  thin  film  heatinq  rate  data  and  illustrates 
very  strikingly  the  "burst"  phenomena  associated  with  trans- 
ition. His  data  very  neatly  illustrate  the  transitional  RMS 


121 


peak  which  also  has  been  observed  in  experiments  usinq  micro- 
phone transducers. 

Examples  of  thermocouple  transducer  assemblies  used  in 
these  applications  are  shown  in  Figure  53;  the  majority  of  such 
devices  which  have  been  used  in  flight  tests  were  designed  and 
fabricated  for  the  test  and  are  not  available  as  off-the-shelf 
devices.  Typically,  the  thermocouples  are  placed  in  a bridge 
circuit  and  connected  to  an  amplifier  which  provides  the  output. 

The  Gardon  gage  in  Figure  53(a)  depends  for  its  output 
on  the  slower  rate  of  heat  transfer  from  the  center  of  a thin 
diaphragm  to  the  edges  which  are  in  contact  with  the  wall  at  a 
lower  temperature.  The  rate  of  response  depends  on  the  diaphragm 
thickness  and  conductivity  and  on  the  center  conductor  diameter. 
This  technique  has  been  especially  successful  in  wind  tunnel 
experiments  wherein  the  pressures  and  temperatures  are  both 
lower  than  in  ICBM  reentry.  It  can  be  seen  that  if  the  foil  is 
thin,  the  device  is  very  fragile  and  easily  damaged. 

The  differential  depth  thermocouple  pair  (or  "AT")  seen 
in  Figure  53(b)  has  produced  excellent  flight  test  results. 

The  General  Electric  Company  has  measured  transition  with 
response  times  shorter  than  10  ms  with  this  approach  which  was 
used  on  PVM  and  STM  instrumented  RV's.  Failure  occurs  during 
reentry  below  transition  altitude  when  the  forward  connection 
is  ablated  away. 

The  final  depiction  in  Figure  53(c)  is  a generalized  view 
of  a calorimeter  which  measures  the  integrated  surface  heat 
flux.  The  calorimeter  output  functions  slowly  when  compared 
with  the  Gardon  or  AT  types  but  is  easily  calibrated. 
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In  addition  to  the  above,  Raman  (Reference  47)  used  a 
very  fast  thin  film  thermocouple  which  responds  essentially  to 
the  surface  temperature.  Mechanical  features  of  the  transducer 
he  used  are  shown  in  Figure  54,  which  was  taken  from  his  report. 

Thin  films  have  been  otherwise  applied  to  measurement  of 
fast-fluctuating  heating  rates  by  using  them  in  resistive-type 
transducers  at  the  surface  of  models.  Owen  (Reference  43,  see 
also  Figure  50  of  this  report)  and  Demetriades  (Reference  45) 
both  have  made  excellent  measurements  with  thin  films  in  wind 
tunnels.  Owens  stated  that  "more  reliable  results  will  always 
be  obtained  using  thin-film  gages"  because  it  is  easier  to 
distinguish  transition  characteristics  from  thin  film  than  it 
is  thermocouple  data. 


Foil  constantan  or  other  metal 


Small  diameter  Cu  wire 


Heat  sink 


Output 


a)  FOIL  OR  GARDON  TYPE 


Cuz  „ 


Thermocouple  output 


b)  DIFFERENTIAL  DEPTH  THERMOCOUPLE  OR  AT  TYPE 


Heatsink  e.g.Cuslug 


Thermocouple  output 


c)  CALORIMETER  TYPE 


FIGURE  53  THERMOCOUPLE  TECHNIQUE 
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FIGURF  54  CRGSS-SECTIONAL  SCHEMATIC  OF  TCS  SERIES  PROBE 
( REFERENCE  16  ) 


4.0  CONCLUSIONS  AND  RECOMMENDATIONS 


4.1  Conclusions 

Analysis  of  the  effects  of  asymmetrical  boundary  layer 
transition  on  conical  reentry  vehicles  has  shown  that  the  most 
important  effect  of  boundary  layer  transition  on  the  vehicle 
forces  and  moments  is  the  boundary  layer  thickening.  This 
effect  causes  an  increase  in  the  normal  pressure  acting  on  the 
vehicle  due  to  the  boundary  layer  displacement  effect.  Using 
a postulated  leeside  favored  transition,  the  induced  pressure 
caused  by  the  boundary  layer  displacement  effect  produces 
incremental  forces  and  moments  acting  opposite  to  the  inviscid 
forces  and  moments.  Since  the  inviscid  forces  and  moments, 
coupled  with  a positive  vehicle  static  margin,  produce  the 
necessary  static  stability  and  angle  of  attack  convergence 
during  reentry,  the  induced  pressure  caused  by  leeside  boundary 
layer  transition  can  be  expected  to  reduce  this  static  stability 
and  may  even  be  a source  of  angle  of  attack  Lvergence  during 
the  boundary  layer  transition  event. 

Comparisons  between  the  wall  shear  stress  and  the  induced 
pressure  caused  by  the  boundary  layer  displacement  effect 
during  transitional  and  fully-developed  turbulent  boundary 
layer  flow  has  shown  that  the  boundary  laycr  displacement 
effect  is  of  significantly  greater  importance  than  wall  shear 
stress.  The  wall  shear  stress  produces  moments  which  tend  to 
decrease  the  static  stability  during  reentry  below  angles  o' 
attack  of  1 . 5 degrees.  Above  1.5  degrees,  wall  shear  stress  is 
found  to  increase  vehicle  static  stability.  Moments  caused  b) 
the  boundary  layer  displacement  effect  are  approximately  a 
factor  of  eight  larger  than  moments  derived  from  wall  shear 
stress  near  zero  degrees  angle  of  attack. 

The  effects  of  surface  mass  transfer  due  to  carbon  phenolic 
ablation  were  found  to  be  relatively  unimportant  with  respect 
to  the  more  significant  boundary  layer  displacement  effect. 
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Surface  mass  transfer  contributes  only  approximately  10  percent 
to  the  net  boundary  layer  displacement  thickness  for  the 
transitional  and  fully-developed  turbulent  boundary  layers. 

If  alterations  in  the  static  stability  of  reentry  vehicles 
are  the  source  of  the  observed  transient  behavior  durinq 
boundary  layer  transition,  then  this  analysis  indicates  that 
small  anqles  of  attack  are  essential  prior  to  the  boundary 
layer  transition  event  before  transient  vehicle  instabilities 
and  angle  of  attack  divergence  can  occur.  This  is  because  the 
incremental  force  and  moment  perturbations  caused  by  asymmetrical 
boundary  layer  transition  are  significant  with  respect  to  the 
inviscid  forces  and  moments  at  angles  of  attack  of  approximately 
1 degree  or  less.  The  larger  inviscid  forces  and  moments 
which  occur  at  the  larger  angles  of  attack  are  sufficient  to 
overcome  the  asymmetrical  boundary  layer  transition  perturbations. 

This  analysis  is  incomplete  in  regard  to  the  accurate 
specification  of  the  dependence  between  the  boundary  layer 
asymmetrical  transition  pattern  and  vehicle  angle  of  attack. 

Many  investigators  have  attempted  to  determine  this  dependence 
experimentally  using  wind  tunnel  simulations.  For  instance, 
several  investigators  have  found  a strong  dependence  between 
the  boundary  layer  transition  asymmetry  pattern  and  both  angle 
of  attack  and  body  bluntness  ratio.  Wind  tunnel  simulations 
of  boundary  layer  transition  effects,  however,  must  be  viewed 
with  some  suspicion  due  to  the  great  disparity  in  the  observed 
transition  unit  Reynolds  numbers  between  flight  and  wind 
tunnel  environments. 

A literature  survey  has  been  completed  with  the  primary 
purpose  to  ascertain  and  assess  methods  which  may  be  applied 
to  measure  pressures  and  determine  location  of  transition  of 
boundary  layers  on  flight  test  reentry  vehicles.  Techniques 
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used  in  actual  fliqht  and  in  wind  tunnel  tests  were  souqht 
and  evaluated  for  the  fliqht  test  application. 

The  conclusions  reached  concerning  available  transducers 
may  be  summarized  as  follows: 

1)  To  measure  vehicle  forces,  static  pressure 
measurements  with  bandwidths  up  to  100  Hz 
desirable.  The  test  results  found  in  the 
survey  of  this  measurement  showed  that  sizable 
errors  are  encountered  in  some  instances  which 
are  related  to  the  interaction  of  the  trans- 
ducer with  the  boundary  layer  flow.  Flush- 
mounted  transducers  and  pitots  inserted  in  the 
flow  to  be  monitored  provided  the  most  accurate 
data  with  errors  less  than  a few  percent. 

Other  installations  which  used  ports  or  inset 
transducers  gave  good  data  during  laminar  flow 
but  produced  erratic  results  during  and  after 
transition.  Errors  of  20  to  30%  were  observed 
in  measurement  of  static  pressure  during 
turbulent  flow  on  flight  vehicles. 

2)  Acoustic  pressure  in  boundary  layers  is  a 
reliable  measurand  of  transition  location. 

Several  experimenters  used  the  noise  from  the 
boundary  layer  and  achieved  good  results,  in 
some  cases  in  the  presence  of  high  level  noise 
backgrounds  in  wind  tunnel  tests.  Both  ported 
and  back-surface  transducers  were  reported; 
back-surface  transducers  have  the  advantaqe  that 
no  penetration  of  the  heatshield  is  required. 
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3)  Temperature/Heating  Rate  Measurements.  Tempera- 
ture measurement  devices  provide  excellent 
transition  location  data  in  moderate  environments 
and  function  reasonably  well  in  the  more  severe 
flight  environment.  Because  of  their  inherent 
heat  integrating  (or  calorimeter)  characteristics, 
rapid  response  is  difficult  to  achieve  especially 
where  measurements  are  to  be  made  of  an  inter- 
mittent turbulent  boundary  iayer.  This  is  due 
mainly  to  heat  diffusive  time  in  heatshield  or 
transducer  materials  being  much  longer  than  the 
changes  in  temperature  at  the  heatshield  surface. 

Experimenters  preferred  thin  film  resistive 
devices  over  thermocouples  if  the  application 
is  permitted.  Flight  test  data  taken  with 
thermocouples  showed  that  transition  beginning 
is  detectable  reliably,  but  later  events  are 
not  detected.  It  should  be  noted,  however,  that 
temperature  is  not  directly  related  to  forces  on 
the  body. 

4.2  Recommendations 

In  order  that  the  effects  of  transition  asymmetries  can  be 
fully  understood  it  is  recommended  that  fast-response  transition 
monitors  be  used  to  fully  map  the  development  of  turbulence  in 
a flight  test  environment.  These  data  would  then  be  compared 
with  the  vehicle  dynamic  response  data  to  determine  the  degree 
of  correlation  between  trajectory  perturbation  and  transition 
asymmetry. 
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In  addition  the  following  recommendations  are  submitted. 


1)  That  new  techniques  be  investigated  for  measure- 
ment of  static  pressure  with  an  accuracy  goal  of 

2 

0.3  n/m  during  the  transition  period.  This  is 
a technically  challenging  task  because  of  the 
reentry  vehicle  boundary  layer  environment,  but  is 
not  beyond  transducer  state-of-the-art.  This 
measurement  need  only  to  have  bandwidths  from  0 
to  100  Hz. 

2)  That  additional  emphasis  be  given  to  boundary 
layer  acoustic  noise  measurements,  especially 
those  which  can  be  made  from  the  back  surface  of 
the  heatshield,  not  contacting  the  Hypersonic 
flow  of  the  boundary  layer.  There  was  some 
optimism  in  the  few  reports  of  this  type  of 
measurement  that  transition  asymmetries  can  be 
readily  detected  with  the  noise  generated  by 
turbulence . 
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